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Abstract 



r- 1 ' Tau functions expressed as fermionic expectation values [14] are shown to provide a 

Qh natural and straightforward description of a number of random processes and statistical 

models involving hard core configurations of identical particles on the integer lattice, like 
a discrete version simple exclusion processes (ASEP), nonintersecting random walkers, 
lattice Coulomb gas models and others, as well as providing a powerful tool for combina- 
torial calculations involving paths between pairs of partitions. We study the decay of the 
initial step function within the discrete ASEP (d-ASEP) model as an example. 

> 

J£} ; 1 Introduction 

■ Random processes and integrable systems are often regarded as opposites extremes in the study 
of dynamical systems, since the first, by its very definition involves indeterminacy and prob- 
abilistic considerations, while the second represents deterministic evolution with the highest 
degree of coherence in time and space. However if we consider the simplest example of motion 
of a free quantum particle we find that this simplest model combines integrability with random- 
ness. The Brownian motion and the motion of quantum particles may be described by similar 
tools, say, in terms of correlation functions and integrals over paths. 

If we consider quantum integrable systems, beginning with the earliest studies on spin sys- 
tems solvable via the Bethe ansatz, it is clear that, like all quantum systems, these involve prob- 
abilistic considerations, involving averaging and statistical correlations, and therefore such sys- 
tems combine integrability with randomness. In fact, many are solvable using the same methods 
as certain related statistical systems since the Hamiltonians may be embedded in the same fam- 
ily of commuting operators as the transfer matrices governing solvable statistical models. Later, 
starting with [21], [22], the success was achieved in applying methods of Bethe ansatz to var- 
ious driven-diffusive and non-equilibrium systems like asymmetric simple exclusion process 
(ASEP). 
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We shall focus on the relations of such systems to classical integrable systems rather than 
to quantum ones. 

It was first shown in the work of McCoy, Tracy and Wu on the Ising model that certain cor- 
relation functions in solvable statistical models satisfy Painleve equations, which naturally arise 
in integrable systems under scaling reductions. In the series of papers by the Kyoto school, a 
systematic approach to integrable systems based on fermionic constructions of tau functions 
was developed, in which remarkable links between the quantum and classical interpretations 
appeared. Classically, a tau function may be seen as a sort of universal "potential" for families 
of integrable equations like the KP and TL hierarchies. (In the Hamiltonian setting, it may be 
interpreted as Hamilton's principle function evaluated along the integral manifolds mapped out 
by a complete set of commuting flows.) At the same time, it has an interpretation as a certain 
correlation function for models of free fermions. Although in the classical integrable systems 
context this formulation at first appeared to simply be an interesting reinterpretation of the 
previous approaches to classical integrable systems, based on isospectral flows of linear oper- 
ators, inverse spectral theory and the Zakharov-Shabat's dressing method [66], it subsequently 
turned out to be a powerful tool linking soliton theory with many other fields of physics and 
mathematics. 

Tau functions have long been known to play a central role in the remarkable links found be- 
tween random matrix models and integrable systems [19], [16], [3]. Subsequent applications of 
tau functions were also found in combinatorics, probability theory, statistical physics and quan- 
tum chaos. (See, e.g., the series of papers by C.Tracy, G. Widom, A.Okounkov, A. Borodin, 
K.Johansson, J. Baik, P. Forrester, P. Zinn- Justin, M. Adler, P. van Moerbeke. In particular 



see [63], [46], [48], [9], [10], [11], [12], [33], [6], [15], [16], [17], [18], [4], [5]). Some of these 



links are quite subtle, while others are still not completely understood. They appear however 
to be quite central to the seemingly endless new applications of ideas and methods originating 
in the modern theory of integrable systems to these other domains. 

In this work, we present some direct links of the fermionic construction of tau functions 
as developed in [14], [32] with certain types of random processes and combinatorial problems 
that arise in relation to partitions. As in [14], [32], we use the language of free fermions, Maya 
diagrams and Young diagrams. This "quantum- like" formulation turns out to be very fruitful in 
providing natural links between many stochastic systems and classical integrability, as well as 
suggesting new methods for approaching the combinatorial aspects common to these problems. 
Suitably interpreted, it will be shown to provide a natural and straightforward description of 
a number of random processes and statistical models. This approach is based on relating a 
natural combinatorially defined basis for the fermionic Fock space to states of the random 
system. To each basis Fock vector is associated, both a. ID configuration of a lattice gas (via 
its "Maya diagram") and a Young diagram, which may be viewed as a region of 2D hard core 
particles. This provides a fermionic setting for well-known models like asymmetric simple 
exclusion processes, ID- lattice log gases subject non-intersecting lattice walkers and models 
of fluctuating interfaces. Other work along similar or related lines may be found, e.g. in the 
work of Okounkov, Olshansky, Borodin, Reshetihin, and others. 

The specific random process or equilibrium statistical model is determined by specification 
of the tau function and the interpretation of time variables. A central tool consists of the well- 
known [60], [61,62] Schur function expansions of tau functions, 
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where the sum is over all pairs of partitions A, [i, which may be deduced very simply from the 
fermionic representation. We show that the tau function r A ijU may be interpreted as a generating 
function for the (unnormalized) transition probabilities between states of certain models of 
stochastic systems. After a suitable normalization, the procedure also leads to discrete analogs 
of models of random orthogonal matrices. 

In the theory of integrable systems the variables t and t are known as "higher times" de- 
scribing the simultaneous evolution under commuting flows, and providing solutions to inte- 
grable "hierarchies", such as the Toda lattice TL and KP systems. In the present work, these 
variables are either chosen to have specific values, thereby defining the statistical weight of 
a given configuration, or just play the role of formal expansion parameters which, through 
the Schur function series allow the tau function to be interpreted as a generating function for 
transition probabilities. 

We use the generating tau function for constructing models of random walks of non-intersecting 
("hard core") particles along one dimensional lattice. We obtain two types of models intro- 
duced by M. Fisher [15]: (a) random turn models and (b) vicious walkers (whose particular 
case M.Fisher called lock step models). In models (a) at each tick of the clock a randomly cho- 
sen walker takes a random step, while in models (b) a certain number of walkers take random 
step in each time instant. In both types of models each site may be occupied by only one walker 
at the same time. 

The alternative approach to the same models is to consider the related Young diagrams, 
which are 2D figures which are in one-to-one correspondence with ID configurations of hard 
core particles. This yields models of oscillating domains. 

Models (a) will be related to the repeating actions of elements of the so-called gl(oo) Lie 
algebra to Fock vectors and tau functions. The number of these actions is identified with the 
discrete time T parameter of the stochastic process. Models (b) will be related to the repeating 
action of elements of the related GL(oo) group. The important role of group GL(oo) and its 
algebra in theory of integrable systems was discovered in [14], [32]. 

The structure of the paper is as follows. 

Second part of the introductory section is devoted to the general conception of our approach. 
Here we introduce a time-dependent normalization function Z v (j) which plays a central role 
in the approach. In section 2 we consider models of type (a). After necessary preliminaries we 
come to subsection 2.3 where we relate a graph to each Toda lattice tau function and consider a 
certain random turn walk along this graph. Here notions of path and of permutation factor of a 
path are introduced. In subsection 2.4 we show that each tau function is a generating functions 
which counts differences of weights for two types of random motion along the related graph, 
which are notable by a value of a permutation factor of their paths. In subsection 2.5 we 
consider known models of Coulomb lattice gas in thermodynamical equilibrium, see [17] for 
a review. The Coulomb potential is modified in a way that the vacuum configuration (the step 
function presenting the Fermi surface) has minimal energy. We point out a two fold relation of 
these models to classical integrable systems. In subsection 2.6 we introduce a model of random 
turn walkers in external potential which may be viewed as a discrete-time version of ASEP (d- 
ASEP) with variable site-depending hopping rates and compare related normalization function 
with the partition function of a model of Coulomb gas in thermodynamical equilibrium. We 
show that Z (t) may be viewed as a certain correlation function in a discrete model of a random 
ensemble of orthogonal T x t matrices. For constant hopping rate we find the asymptotic 
configuration in the large T limit. We found that the limiting shape coincides with one obtained 
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by Vershik and Kerov in [38] in the study of a limiting shape of random partitions with fixed 
large area which are distributed according to the so-called Plancherel measure. In our case the 
size of asymptotic configuration depends on the hopping rate as it is given by formula (92). In 
subsection 2.7 we introduce additional dependence of the hopping rates which now depend on 
a mutual position of particles via the modified Coulomb potential with an arbitrary charge and 
show how it corrects the size of the asymptotic configuration. In section 3 we consider models 
(b). In subsection 3.1 we show how the so-called determinantal ensembles [13] appear as a 
result of actions of elements of GL(oo) group on tau functions. In subsection 3.2 we mark the 
convenience of usage of multi-component fermions which allows in a simple way to construct 
fermionic representations for models of vicious walkers in an external potential and via Wick's 
theorem yields answers for correlation functions in form of determinants. In subsection 3.3 
we consider models of "random layering" and models of "rocks" (figures which generalize 
diagrams of the so-called 3-D partitions) which we believe were basically known in certain 
versions, for instance, as the so-called polynuclear growth model. Here again we introduce 
the external potential which results in a variable site-dependent rate (this is the hopping rate in 
ID picture and a gluing rate in 2D one). The normalization function of some of these models 
coincides with the tau function of Pfaff lattice [4], [35]. 

Except the consideration in subsection 2.6 we restrict ourselves to the presentation of var- 
ious models of random processes arising from integrable systems more than to the analysis 
of random models themselves which is mainly an analysis of the large time behavior of dif- 
ferent correlation functions. As it was very successfully demonstrated in series of papers by 
A.Okounkov the usage of fermionic language may be quite helpful for such purposes. This part 
of work will be published separately. 

1.1 Action on Fock vectors. Time-dependent normalization function Z v (t) 

We refer the reader to Appendices A.l- A.3 for a compendium useful facts and conventions 
regarding fermions, partitions and Schur functions that are used repeatedly in what follows. 
For present purposes, this involves only charged, one-dimensional free fermions satisfying the 
standard canonical anti-commutation relations 

1,1) ■ J)L <\r fifj + f S fi = 0, fifj + fjfi = 0, ijez (2) 

These act on a fermionic Fock space whose "vacuum state" |0) is, like the Dirac sea, annihilated 
by all negative component creation operators f { and positive component annihilation operators 

fi 

/_ i -i|0>=/ i |0) = 0, *>0 (3) 
More generally, we have the "charge n" vacuum states \n), 

\n) = fn-i ■ ■ ■ /o|0) for n > 0, \n) = f n ... /_i|0> for n < 0, (4) 

where the integer n denotes the Dirac "sea level". We may visualize the level - n Dirac sea as 
a diagram consisting of a set of integer sites on a vertical axis, with all sites below n occupied, 
while the n-th site and all those above it are empty, as in Fig. 1 . 

A complete set of basis vectors for the Fock space may be associated with Maya diagrams, 
which we shall denote by Greek letters A, v, . . . , to be more precise, by two letters, Greek and 
Latin, the first denotes the configuration of particles, the second denotes the level of the Dirac 
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1 . Dirac sea of level n 2. One particle one hole 

Figure 1 : Maya diagrams 

sea, which we perturb by placing a finite number of particles of the Dirac sea to the upper 
positions. In this way all basis vectors of Fock space may be obtained. Since no two particles 
can occupy the same site, each Fock basis vector may be encoded by a partition, these are the 
ordered sets of weakly decreasing numbers, say, v — (v\, v 2 , . . . ) where v\ > v 2 , ■ ■ ■ > 0, and 
by a level, say n, related to the coordinates of the particles, hi > h 2 > ■ ■ • , as follows: 

V{ = hi + i — n (5) 

We number (identical) particles in downward direction, the highest one be the first. Provided 
we view each Fock vector as being created by a perturbation of vacuum ('initial') vector, one 
may say, that i/j measures the deviation of a current ('Lagrangian') coordinate, hi, of the i-th 
particle from its initial value, n — i, in the vacuum state. 

The partitions are very useful and conventional notations, that is why we need them. There 
are two suitable way to figure partitions: via Maya diagram, which is just configuration of the 
particles h = (hi, h 2 , ■ ■ ■ ) placed on the vertical line, and via Young diagram which is a set 
of v = (i/i, v 2 , . . . ) related to h = (hi, h 2 , . . . ) by (5). They also yields links with the topic 
of growing interfaces, since the diagram of each partition (Young diagram) may be viewed as 
a discrete model of a two-dimensional region. A set of partitions placed one above another 
in three-dimensional space may be a discrete model of 3D region, and so on. A particular 
case of such 3-D figure, called plane partition, was considered in the papers of Okounkov and 
Reshetikhin [48], [49] as a certain statistical model. 

In the present paper we consider a sequence of Fock vectors obtained by application of 
operators d 1 ^ 1 to an initial basis Fock vector, say, \u),as follows 

|i/,n)->o 1>0 |i/,n)->o 2 'V >0 |i/,n) -> > o 1 ^ 1 ■ ■ ■ o lfi \u, n) (6) 

This sequence may be viewed as a T-step process on the space of configurations of hard 
core particles (ID lattice gas) where each configuration is related to some basis Fock vector 
and may be visualized via the so-called Maya diagram, see the picture. The number T will be 
related as discrete time of the process. Alternatively, (6) may be viewed as a T-step process on 
the Young diagrams (a type of 2D regions). 

In what follows everywhere, where it will not be not confusing, we shall omit the label 
related to the level of configuration (or, the same, the level of Dirac sea), keeping only Greek 
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letters, say, A, fj, ,v for configurations of the hard core particles, having in mind the pairs (A, I), 
(/!, m) n). 

The result of the application of 

o(t) := o 1 ' 1 " 1 • • • o 1 ' (7) 

to the initial state v (corresponding to an initial configuration of hard core particles), is a Fock 
vector which is a linear combination of basis Fock vectors (=configurations of hard core par- 
ticles). One of the objects of an interest is the relative weight of a given configuration, say, 
A, 

W^a(t) = (A|o(t» (8) 
with respect to the weight of all possible configurations, 

p ( T ) - jWlj m\ 

p ^ a(t) - (9) 

where the normalization function in the denominator is 

Z„(t)=£<A| (t)|i/> (10) 
A 

the sum ranges over all possible configurations. 
Matrix element 

^i~ 1 = <M° , " J "~V"»> (ID 

will be referred as (one-step) transition weight between a configuration m) and a configura- 
tion | A, /) at moment j. 

If we want to keep in mind the probabilistic interpretation of the process we ask transition 
weights to be positive. 

The completeness of the basis Fock vectors implies ^ m m) (//, m\ — 1. We have 

(XJ\o^...o^\u,n)= J2 oj;ti)-^S V < 12 > 

{juW,mi} 

where one sums over all possible intermediate configurations ^ % \rrS l \ % — 1, . . . , T — 1. 
Each monomial non-vanishing term in the sum is related to a certain set of configurations, 

v,n^ii {1 \m^ -> >■ /x (T_1) ,m (T_1) — > A, Z , (13) 

which will be refereed as a from the configuration u, n to the configuration A, I, having 
length (or, the same, duration) T. The value of the monomial is referred as the weight of the 
path, which is equal to a product of all one-step transition weights (11) along the path. 

Let us notice that it may have sense to evaluate Z„(t) also for the case where the positivity 
condition, 

«r > o. (1 4 ) 

is violated and therefore we have no a probabilistic interpretation of (12) as a random process 
on the space of configurations of the hard core particles. 
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In certain examples considered at the present paper transition weights (11) take the form 

(A, Z|o™- V, m) = e-^M+u^mU-i) ( 1 5) 

The case where the weight of each path depends only on the end points of the path, i.e. on 
initial and on final configurations shall be referred as potential weights. In what follows we 
shall consider only potential transition weights. 

To compare with, let us assign to each configuration of hard-core particles, A, an energy U\. 
What people are interested in statistical physics is the partition function of the system 

£ e~ U > (16) 

A 

(everywhere we shall put temperature times the Boltzmann constant to be one). 

In case of random process the normalization function (10) - which yields the sum of weights 
of all pathes of durations T which start at the configuration v - plays the role similar to the role 
of partition function in thermodynamics. For potential transition weights the normalization 
function takes the form 

Z„(T) = e-^WW) N V (T) (17) 

A 

where N a , v (t) counts the number of pathes of duration T. Given T, this factor may be also 
included into the exponent to define a partition function of a T-dependent statistical model. 

At the present paper we shall focus to figure out the links of random processes (6) to clas- 
sical integrable systems, random matrices and growing surfaces. Benefits of these links will be 
presented in our forthcoming papers [29], [30]. 

Our main example is a discrete version of asymmetric simple exclusion process (d-ASEP), 
where we present the answer for the decay of the step function. In short we shall consider 
random walk of N non-intersecting particles and other examples. 

The mentioned links of random processes with classical integrable systems are two-fold. 
The first one is the usage of series (1) as generating functions for transition weights between 
different configurations of the hard-core particles. Then higher times t and t are just formal 
parameters. There are different tau functions, related to A, B, D and C root systems; the way 
of construction of these may involve different type of fermionic (or, bosonic) operators. At the 
present paper we deal with A type, which may be realized in terms of fermions (2). We show 
that A-type tau function may generate certain classes of transition weights (11). We consider 
two ways of generating weights, see below. This is one link between random processes and 
integrable systems. 

The second link is the fact that the normalization functions Z (t) may be related to a tau 
function of some other classical integrable systems. It occurs if we consider deformations of 
transition weights with a set of deformation parameters t — (ti, t 2 , ■ ■ . ) 

(X,l\o j,j ~ 1 \fj 1 ,m) -> (X,l\o i ' j - 1 (t)\fx,m) = e -^.«CJ,*)+^,m(j-i,*)^, m) (18) 

specified later by (58) and caused by transforms of operators o 7 '- 7-1 of form 

ojj-i ^ ow-i(f) = e ^o(t) jJ-i e -Ho(t) > j = 1 , . . . , t (19) 
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(operators H are defined in Appendix A. 3). These deformation parameters may be recognized 
as higher times of different underlying integrable equations. These are dual integrable hierar- 
chies in sense of papers [50] and [54] but not exactly 4 . In many cases, under flows (19), the 
normalization function Z (t, t) of (10) is a tau function of B-type, namely, the B type in the 
realization of Kac-van de Leur, see [34] (also known as Pfaff lattice independently introduced 
as a tool for study of random matrices see [4], [40]). Thus, what we mean is a sort of duality 
between A and B type of equations. 

The last remark is that we may consider a sort of coupling of random processes with statis- 
tical models on the space of configuration, 

Z = J2 e ~ q2Ex > (20) 

A 

where q is some constant (charge) as follows. We use the energy g 2 E A of the model to modify 
each transition weight: 

(Al^-V) - e-« 2E ^ 2E " (A|^-V>, J = 1,2,...,T (21) 

Then 

Z u (T) = J2(Mo(T)\v) -> Yl e~^ W *» <AKT)H (22) 

A A 

As we shall show in subsection 2.7 this procedure may keep relations with integrable sys- 
tems if q 2 E\ is chosen in a special way. 

2 Action on Fock vectors via gl(oo) algebra 

Maya diagrams 

Let us consider an infinite set of vertexes labelled by integers. Each vertex will be sketched 
as either a white or black disk. Following Sato, we shall call it Maya diagram, (after a Maya 
game in which small stones are moved to occupy free sites according to a certain rule). In what 
follows we shall describe certain rules for motion of "stones"- hard core particles. The Maya 
diagrams for the vacuum state are presented in Figure 1.2. 

We may consider each black ball as a particle and each white ball as the lack of a particle 
at this site (empty site). 

The Maya diagram is a way of visualization of basis Fock vectors for fermions on the circle, 
sites correspond to Fourier components of the fermions, no two fermions can occupy the same 
site; therefore, there is one site for at most one particle (these are what is called hard core 
particles). Sites with positive large enough numbers are supposed to be free, while sites with 
large enough negative numbers are supposed to be occupied, that is the related Maya diagrams 
possesses this property: all sites to the up of a certain number are all white, while sites to the 
down of a certain number are all black. The Maya diagram where all sites below n-th one are 
occupied, and all sites above n — 1 are free, is called the vacuum Maya diagram of level n, see 
the Figure 1 . It represent the so-called Dirac sea of level n. 

4 This duality might be related in a general notion of duality of integrable systems developed in [2], [24], [25], 

[7] 
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Each Maya diagram may be obtained from the vacuum diagram of a certain level, say, n, 
via placing of certain number, say, k, of black balls (particles) to higher empty sites. We say 
that each Maya diagram has level n. 

Basis Fock vectors, partitions and Maya diagrams 

By its origin Maya diagrams is a visualization of vectors in the fermionic Fock space, see 
Appendix A.l. Each black ball (a particle) located at a site % indicates that the site is occupied. 
Let us introduce the following basis vectors in the fermionic Fock space 

\X,Tl) = (-1) E ' =1 ft /n+a fc /n-l-/3 fc " ' fn+aJn-l-pM) i a ± > ■■■>a k >0, ft > • • • > ft > 

(23) 

where A = (a±, . . . , afc|ft, ■ ■ ■ , ft) is the Frobenius notation for the partition A, and n is the 
level of Dirac sea, see Introduction and Appendix A.l for more details. 

This vector is one-to-one related to the Maya diagram of level n, where particle are situated 
as follows. A particle number % (black ball) counted from the top has a coordinate n + a i5 
% — 1, . . . , k. Thus, the upmost particle has the coordinate a\. Underneath of the sea level all 
sites are occupied except free sites (holes) which we count from the bottom. Hole number i has 
a coordinate n — 1 — ft, thus the down-most hole has the coordinate n — 1 — ft. 

Notice that ft + 1 is the length of partition lambda, 

One may give an alternative description of coordinates of particles on the Maya diagram 
related to the vector |A, n). Take any integer N which is not less than £(X). Then all site with 
coordinates n — N — 1, n — N — 1, . . . are occupied. On the sites n — N, n — N + 1, . . . there 
are N particles which occupy sites with coordinates 

hi = \i-i + N, i = l,...,N, N>£(\) (24) 

All other sites including and above n — N are empty. 
Dual vector is 

(\n\ := Afcl/.^/^ . . . f n ^J n+ctk ■ ■ ■ (25) 

One can verify the orthonormality condition 

(\,l\v,n) = 5i, n 5\,v (26) 

Along these notations we shall identify partitions, Maya diagrams and basis Fock vectors. 
Being identical to basis Fock vectors, Maya diagrams may be multiplied by numbers and may 
form (formal) linear combinations. 

Action of gl{oo) on Maya diagrams 

It was old observation that quadratic operators of type fifk, i, k = 1, . . . , N, form gl(N) Lie 
algebra. In the papers of Kyoto school it was studied within symmetry analysis of integrable 
equations like Toda lattice equation. One of the observations made in the wonderful series 
devoted to hierarchies of integrable equations (in particular, see [14], [32]) is that N = oo case 
needs a special consideration. When there are an infinite number of fermionic modes one needs 
normal ordering. Then, natural combination are 

E i>k := fif k - (0\fj k \0) - oo < i, k < +oo (27) 

which form an algebra called gl (oo) (the algebra of infinite matrices with central extension). 
We do not need properties of this algebra. Let us only note that such algebras play essential 
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role in the construction of hierarchies of integrable classical (not quantum) equations, where the 
algebra of infinite matrices appears as an underlying symmetry algebra of integrable equations. 
This algebra naturally acts on the fermionic Fock space. 

For our purposes we are interested in the action of E ijk on the basis Fock vectors, which 
may be viewed as an action on Maya diagrams. The point is that this action transforms one 
Maya diagram to the other one multiplied by a sign factor. Then any product of operators (27) 
do the same. 

Keeping in mind definition of basis Fock vector and applying (27) one can verify that the 
action of E ijk on Maya diagram is as follows. 

First, let us notice that the second term in the definition of operator E ijk is vanishing for % ^ 
k. Let us consider this case first. Each E ijk acts on a Maya diagram trivially (as multiplication 
by zero of the related basis Fock vector) in two cases. They occur either if k-th site is empty, 
or/and if i-th site is already occupied by a particle as on Figure 2.1. Otherwise E iyk action 
places black ball from its A;-th site to the % site, multiplying the new Maya diagram by (— l) Cik , 
where c ik is a number of black balls between i-th and /c-th sites, in other words, the number of 
particles which the particle originally located at A;-th site "jumps over" when moving to its new, 
i-th, site, see Figure 2.2. An arrow on this figure shows that the particle hope from the site k to 
the site i. 



E; 



i,k- 



• i 



= 0, 



E; 



i,k" 



o k 



= 0, 



o 1 i 

9ik 



-i.k- • = (-D 
• k 



v. 



o k 



1 . Elimination of Maya diagrams 



2. Nontrivial action 



Figure 2: Elementary gl(oo) actions on Maya diagrams 



The action of operators E i:i on Maya diagrams may be called testifying action. They acts 
trivially on the vacuum Maya diagram of level 0. For i < the operator E iti eliminates any 
Maya diagram which (similar to the vacuum diagram) contains a particle at the site % < 0. If 
the site % is empty then E i:i acts as the multiplication by unity. For % > the operator E i:i 
eliminates any Maya diagram which (as the vacuum diagram of the level 0) contains the white 
ball at the site i. If there is a particle at the site i then E iyi acts as multiplication by unity. 

Because Maya diagrams are in one-to-one correspondence with basis Fock vectors we can 
act on them by sums and products of operators (27). 

The element 

plays a special role. One can easily check that operator Q acts a Maya diagram of a level n by 
the multiplication by n. It is called the charge operator. 
Graphs 
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Let us take vertices of a Maya diagram as vertices of a graph, where vertices are labelled 
by integers and each pair of vertices may be connected by a pair of opposite directed arrows. 
Given 

A = J2^,kEi, k egl(oo), (29) 

we assign the weight a iyk to the arrow which starts at the vertex k and ends at the vertex i. For 
the sake of simplicity, we shall not consider diagonal terms, E iti . Also let us not draw arrows 
of zero weight. In this way we obtain some weighted graph, related to given A =G gl(oo). For 
instance, for 



A= ( E i+hi + E i,i+i) (3°) 



=n-2 



we obtain the graph on Figure 3.1. 
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1. Graph 2. 3. 4. 5. 

Figure 3: Maya diagrams projected on graphs 

In the present paper we are interested in motion of particles along arrows of the graph 
whose vertices coincide with vertices of Maya diagram. Having in mind this picture one may 
forget about all vertices of Maya diagrams which are not linked by vertices. For instance, the 
Maya diagrams depicted on Figures 3.2 and 3.3, projected on the graph depicted on Figure 3.1, 
respectively, are as on Figures 3.4 and 3.5. If the number of arrows is finite we shall ignore the 
rest part of a Maya diagram. 



2.1 Young diagram versus Maya diagram. Interface between 2D area 
versus configuration of particles on the ID lattice 

For discrete systems the difference between ID and 2D systems of hard core particles is 
rather conventional in the following sense: any ID lattice gas configuration may be converted 
into a Young diagram viewed as an example of a condensate of hard core 2D particles. 

Let us figure out an example: 

These two figures show the relation between the motion of particles on Maya diagrams and 
a motion of the interface. 

One may say, that each hop of a ID particle upward is related to an increasing of the area of 
the corresponding Young diagram (a condensation of 2D particles to the Young diagram), while 
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1 . Upward/downward hop 2. Adding/removing a box to a Young diagram 

of a particle on a Maya diagram 

Figure 4: Equivalent representations of elementary steps 

a hop downward results in the inverse process of 'evaporation of particles' and the decreasing 
of the area the Young diagram. (A figure similar to Fig. 4.2 may be found in [14] without 
relations to the topic of the present paper: to stochastic motion.) 

2.2 Generating vector for Maya diagrams 

There are vectors, coherent states, (n|e H ( t ) and e H ^\n), which depend on parameters t = 
(ti, t 2 , . . . ), see Appendix A. 3, and which generate all Fock space with fixed given sea level n, 

{nle H(t) = J2(X,n\s x (t), e S ^\n) = £ \X,n)s x (t) (31) 

AGP AeP 

where S\(t) is the Schur polynomial [42]. The Schur polynomial is a polynomial in many 
variables t±, t 2 , ■ ■ ■ , which are labelled by partitions. The appearance of the Schur polynomial 
is typical for integrable systems. 

2.3 Random turn motion on graphs. Permutation factor and weight of 
path 

Consider the graph related to 

A = a i,kE itk , E itk = fif k , i^k (32) 

which is a set of vertices, some of which may be linked by arrows: an arrow (i, k) starts on a 
vertex k and ends on a vertex i. We assign a weight a i)k to each arrow (i, k). In case a i>k = 
there is no arrow. The matrix a with entrances a i>k is called connecting matrix for a graph. For 
simplicity of consideration we omit diagonal terms. 

We shall consider a discrete dynamical motion along a Maya diagram with an initial dia- 
gram v given by 

\v) -> A\v) -> ► A T \v) (33) 

Here T is viewed as discrete time variable. 
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As we have explained in the Introduction we have different paths connecting two different 
configurations, see (13); each path is related to a sequence of intermediate partitions. 

Consider A T which is a sum of monomials, each one is a product of T terms a ia yka E ia ika , a = 
1, . . . , T. Each path of process (33) is generated by a unique monomial applied to a initial Maya 
diagram. The number T may be referred as duration of path. It may be described as set of con- 
sequent events, single hops. Each time instant, say, a, a particle hops along an arrow (i a , k a ): 
one hop at each time instant. Say, monomial a iT> k T E iT:kT ■ ■ ■ a iljkl E iljkl , provided it acts non- 
trivially on v, at first time instant describes the hop of a particle along arrow k\) and at last 
time instant, T, describes the hop of particle along an arrow (i T , k T ). The product weights of T 
arrows which were crossed by particles, Yla=i a i a ,k a > is called the weight of path. 

Each pair of configurations may be connected by a number of paths each of which has its 
own weight. 

The fact that A T is a sum of weighted monomials means randomness of the motion de- 
scribed by (33). 

Let us call it random turn motion along a graph, given by A. 

The random motion where in each time instant only one chosen at random particle hops is 
called random turn motion. 

To each path we shall assign permutation factor as follows. 

Given acting nontrivially monomial, we can follow each particle involved into the motion 
along graph. Let us enumerate particles located on the initial Maya diagram the from the top 
of diagram, and keep these (personal) numbers of particles on the target Maya diagram (where 
each particle 'remembers' it's given number). Let us re-enumerate particles of the target Maya 
diagram from the top: this may be also obtained by a certain permutation of personal numbers 
which we assigned to the particles in the beginning of the motion. If the sign of the permutation 
is even then the permutational factor of the path is equal to + 1 . If the sign is odd the permutation 
factor of the path is equal to —1. 

Thus, to a each path we assign a weight, a permutation factor, and a duration. 

2.4 Graphs and tau functions. Tau function as counting function 

Consider the random turn motion described in the previous subsection. 

Consider all paths of duration T starting with a configuration described by a partition A' and 
ending on a configuration A. The sum of weights of all such paths with positive permutation 
factor we denote by 



Let us note that permutation sign appears due to the fermionic approach to the problem. 
In problems of random motion of hard-core particles along graphs which know nothing about 
fermions, one may wonder about a total weight, which is 




(34) 




(35) 



W^ x (t)=W^I x (t)+w[tI x (t) 



(36) 



From the previous consideration we have 



(X,n\A T \X',n} = W^ x (t) -W^ x (t) 



(37) 
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By (31) we relate these numbers to the so-called TL tau function [14], [64] in the following 
way 

OO j 

r A (n,t,t) := (n\e H Ve* A e S V\n) = £ £ ^ { W ™xW ~ W$~_> a (t)) s x (t)s x ,(i) (38) 

T=0 A, A' 

where A is 

\ >>.,/,/:, (39) 

i¥=3 

and z, t = (ti, t 2 , ■ ■ ■ ), t = (?!, i 2 , ■ ■ ■ ) are formal parameters, and f i: fa are free fermions. 

oo oo oo 

H(t) = ^2 Hmtm, H(t) = ^2 H- m t m , H m = ^2 fifi+m (40) 
m=l m=l i=—oo 

Thus, we say, that the power of z counts the number of steps. T-step process is 

\v)^A\v)^ >A T \v) (41) 

The variable T is treated as the (discrete) time of the random process described by A. 
The most simple case is random processes where 

W ( x 7_l x (T) = (42) 

for all pairs A' A. Then tau function may be interpreted as counting function for weights 
W / a'^a(t). Then it follows that Wa'^a(t) is a subject to certain discrete bilinear equations 
(Hirota equations), which we shall write down not here. 

However we can weaken condition (42) keeping the interpretation of tau function as count- 
ing function. Consider two examples where particles may hop only to the nearest neighboring 
sties. One example is presented on the figures 5. Figures 5.1 and 5.3 depict graphs where ini- 
tial step function, v — 0, configuration will decay. In figure 5.1 at each time instant a particle 
chosen at random will hop upward along an arrow if the neighboring place is not occupied. 
In figure 5.2, motion is impossible: in this configuration particles are locked. Each finite per- 
turbation of v — will come to v = in a final number of steps. In Figure 5.3 not only 
particles are chosen at random but also arrows are to be chosen at random: if the neighboring 
vertices are free, chosen at random particle may hop either upward or downward. It is called 
non-intersecting random turn walkers, (on this topic see [6], [17], [8] and [3]). These examples 
of random motion keep the condition (42) because it is ordered motion and the signs of each 
path is positive. For the choice of A as in Figures 5.1-5.3, namely, A — Hi, Hi + 
the weight of any path is either zero, or one. It means that tau function is a generating function 
for number of ways to get configuration A starting with a configuration A' in T steps. 

Second example is a walk on a ring, where particles are allowed to hop to the nearest 
neighboring site. How to obtain such graphs is depicted on figure 6: figures 6.3 and 6.4 are 
obtained respectively from figure 6. 1 and 6.2. On figures 6.2 (or, the same 6.4) condition (42) is 
violated because when a particle hops along the arrow linking vertices n and n + 5 permutation 
of particles is odd in case the total number of particles on the ring is even. Thus each path 
containing odd number of hops along arrows linking sites n and n + 5 has negative permutation 
factor. Nevertheless if we assign positive weight to all arrows except arrows linking n and n + 5 
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Figure 5: Random turn motion along graphs 
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Figure 6: For the case depicted in 2 and 4, in order to keep the probabilistic interpretation, one 
has to assign negative weights to the arrows that link n + 5 and n 



sites, tau function will count the transition weights between any pair of configurations. Say, if 
arrows linking n and n + 5 have weight —1, while all other arrows have weight 1 tau function 
counts number of ways linking any pair of configuration on the ring. 

Processes like those depicted in Figures 5 and 6 are also known as discrete versions of 
simple exclusion processes (totally asymmetric (d-TASEP) for figures 5.1, 5.2, and symmetric 
for figure 5.3), which may be considered either on the line as on figures 1-3, or on a ring as 
in the Figure 6. This type of models was introduced by M.E.Fisher. It may be called discrete 
time asymmetric simple exclusion process (d-ASEP) in case the probability rate for hops in 
the opposite directions are different and does not depend upon site numbers. The solvability of 
ASEP was proven by H. Spohn in [21, 22] by method of Bethe anzats and a huge literature is 
devoted to the study of different aspects of this model which, in non-equilibrium physics, plays 
a role comparable to the role of Ising model in statistical physics. 

The problem of evaluation of correlation functions, asymptotic behavior we hope to con- 
sider separately using fermionic approach. 

The dependence of transition probabilities on parameters (weights of arrows) will be con- 
sidered in subsection 2.6. Before this consideration of non-equilibrium system of particles we 
find that a consideration of fermions in thermal equilibrium will be helpful. 
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2.5 Fermions in thermal equilibrium 



Consider the following statistical ensemble: there is an infinite set of levels labelled by inte- 
gers. Each level % has energy Ui and may be occupied by a fermion. In thermodynamics the 
probability to occupy a site % is proportional to exp(— Ui), where U is energy level numbered 
by i (we put the Boltzmann constant times temperature to be one). In case U is a monotonous 
function of the site number we can use Maya diagrams of some level, say n. At zero tempera- 
ture all lower site are occupied up to this level (in solid state physics called Fermi level). The 
bottom limit of Maya diagram is fully packed by particles, while top limit is free of particles. 
A fragment of the picture is sketched out in the following figure 

• i+1 exp(-U i+1 ) 
o i exp (-U-) 

• i-1 expfU^p 

• i-2 exp (-U-_2) 

Fig.l. Fermionic levels where some sites are occupied 



Each configuration of fermions distributed among levels has its weight given by the product 
of the Gibbs factors related to the occupied sites. The contribution of the fragment depicted on 
the figure is e - u i-2-Ui-i-u i+1 ^ 

These are non-interacting fermions. The equilibrium ensemble is described by its partition 
function. 

First, let us introduce the notation which we shall use throughout the paper 

oo 

Ux(n) = ( U ^-i+n - U-i+n) (43) 
i=l 

Notice that for the vacuum configuration of any level n we have U (n) = 0. 
The normalized partition function is then a sum over partitions 

e -Fo = ^ e -U x {n) (44) 
A 

where F is the free energy of the system of non-interacting fermions living on a Maya diagram 
of a level n. 

Remark. Let notice, that we can cut Maya diagram from the bottom via the following 
procedure. We send all energies of sites located bellow, say, site N, to minus infinity. Without 
the loss of generality we take N = 0. Then, we have in particular, 

e U-!-U = o (45) 

In this way we restrict our consideration by finite number of fermions, n. It is available in 
case our system is related to one-dimensional solid state system where there are finite number 
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of particles. In such cases n will be the number of particles equal to Fermi level. If only n 
fermions are involved we have 

n 

) , (46) 

i=i 

where A« — i is the coordinate of the particle numbered by i (i — 1, . . . , n) on a Maya diagram 
of a level n which is cut at zero site. We count particles from the top and their coordinates from 
the bottom. 

We shall consider both cases, N = and iV — > — oo, the first will be described with the 
help of semi-infinite Maya diagram and the second by usual Maya diagram. 

Modified Coulomb potential E A . Now, suppose that fermions pair-wise interact, that is, 
now, the energy of the system is U\(n) plus the energy of pair- wise interaction. We choose this 
interaction as Coulomb interaction (if sites are interpreted as level under Fermi surface then, 
this is Coulomb interaction in the momentum space) 

g 2 log(*-j), (47) 

which describes the repulsion of fermions of a charge q at sites i and j. The problem that 
Coulomb energy goes to +oo for semi-infinite lattice. In order to avoid the divergency of the 
energy for semi-infinite lattice we shall modify it taking the energy of particles in a configura- 
tion A in form 

E A = -log^too) := log n ^f* tH = Xi-i + L, i = l,...,L, L > £(X) (48) 

rii=i hi\ 

Notation s^t^) is taken from the Appendix A.l, see (177). 

Remark 2.1. This expression contains pair-wise Coulomb interaction (47) in the enumerator and an 
external electric potential in the denominator. As one may check this expression vanishes on the vacuum 
configuration for any value of L > 0. Let us mark that L is not a parameter: the energy does not depend 
on the choice of L provided it exceeds the length of a partition. It is reasonable to recall the meaning 
of numbers i, Aj and h\. The number i numerate particles counted in the downward direction. The 
number \ measure the shift of the particle i in the upward direction from it's homesite in the vacuum 
configuration. Given L the number hi is the coordinate of the particle numbered by i which is counted 
from the origin located on the site n-Lon the Maya diagram of level n. The point that the interaction 
(48) does not depend on the choice of this origin which one can send to — oo. It is perfectly adopted for 
the lattice gas with fully packed "bottom" limit. Moreover, one can show that the vacuum configuration 
is the configuration with minimal energy equal to is zero. 

The related free energy will be denoted by F q . 
In this case the partition function of the system is 

e -F B (n) = ^e-tWE^iAiiogCMO = Y j e- u ^ ( SA (too)r 2 (Mi) |A| (49) 

A A 

where we introduce dependence on axillary parameters titi for further convenience. 

There are special cases, q 2 = 1 and q 2 = 2, where (49) may be identified with tau functions 
of quite different integrable systems. (The case q 2 = 4 is also related to integrable systems, 
however we will not need it at present paper). We shall consider these cases separately. 
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For q 2 = 2 the slightly modified series (49) is equal to 

e -MnMM) ;= J2e- U ^ E A+|A|iog(tifi) = J] e -^(f lf -)|A| ( SA ( too )) 2 

A A 

which is a Toda lattice (TL) tau function. (Here c n is a normalization constant chosen in a way 
that the right hand side is equal to unity for t\ — t\ — 0, see (211)). The equality directly 
follows from (31) and (26) if we take into account (210) and finely use (173) and (177). This 
formula is true for both cases N = — oo and N = 0. To obtain the last case we need keep in 
mind relation (45). 

Let us refer to the paper [41] where the Coulomb gas on the lattice was related to the tau 
functions of different type. 

In general, the following vacuum expectation value which is known to be an example of the 
Toda lattice tau function, namely 

r(n,t,U,t) = {n\e H{t) e^ UihIi -^<° UJih e m \n) = c n ^ e~ u *s x (t)s x (i) (51) 

A 

yields the partition function for the system of fermions in a configuration | A, n) whose pair- wise 
interaction (47) is replaced by interaction of all particles via the potential given by 

logs A (t)+logs A (t) (52) 

which, in general, is not of a pair- wise type. The second equality of (51) is derived from 
relations (31),(26) and (210), exactly in the same way as previous formula (50) but without 
specification of Schur functions. The pair-wise type interactions one may obtain via specifi- 
cation of parameters t and t which enter the tau function (51), via formulas (179)-(176) of 
Appendix A. 2. (Such specifications were considered in [52], [53], [54] and [50] in quite differ- 
ent contexts). 

Remark 2.2. To complete the talk about links between the Coulomb lattice gas and tau functions we 
note, that if we choose a special parametrization of the set {Ui,i G Z}: via sets of variables t := 
(ii, ?2j • • • ) and t* := (t_i, i_2, . . . ) as follows: 

oo 

Ui=J2 im tm, i e Z (53) 

m^0 

we obtain (see [50]) that e _p2 is a TL tau function in variables of a certain dual hierarchy. 

Similar parametrization of the parameters U which enter tau function (50) was used in the context 
of completely different problems in [47] and in [1]. In [1] relation (53) has a meaning of the dispersion 
law for fermions in solid state physics, where i plays the role of momentum and Ui plays the role of 
energy of fermions. Authors thanks P. Wiegmann for explaining [1] before it was published. 

The case q 2 = 1 is most important in view of its application in the next section, and is less 
trivial from the point of view of what it may be in integrable hierarchies. Given level of vacuum 
configuration n consider 

6n {U) := (\n\e^>° UifiIi -^<° Uififi e H - 1 \0,n) (54) 

A 

C(X)<n 
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11 e -u hi +u- l+n .™ 



= Cn j2 e-^s^) = Cn j2 E II — )n — - h i) 

A fc=l hi> ->h„>0 i=l l ' i<j 

e(X)<n ~ 

where hi = Aj — i + n. The second equality is obtained from formulas (210) and (206) and 
the last equality from (177) of Appendix A. 2, and where c n is given by (211) (We shall omit 
this constant c n below). Then replacing summation over the cones hi > ■ ■ ■ > h n > by sum- 
mation over non-ordered non-negative hi, . . . , h n which may be done due to the permutational 
symmetry of the terms (and getting the factor (n!) -1 ) we arrive at the partition function for n 
fermions with the electric charge 1: 

1 n —Uh +U- ■+ n 

^ U ) = n\ £ I T ;,,"'" lI^-^- e " Fl(n) (55) 



hi,...,h n >0 i=l i<j 

Notice that we do not need condition (45) because the reduction to the n particle partition 
function is achieved by imposing the condition £(X) < n. 

Remark 2.3. The interesting fact is the following. Let us write the grand partition function 

oo 

e {^U)=Y J ^Qn{U) (56) 

n=0 

where /j, is a chemical potential. Then one can show that this is an example of the infinite-soliton 
tau function of a version of the BKP hierarchy, suggested by V.Kac and J. van de Leur [35] (called 
them "charged" BKP hierarchy). This BKP hierarchy is different from the BKP hierarchy presented 
in [14], [32]. 

The fermionic representation is a specification of the fermionic representation found in [40] in the 
context of the study of Pfaff lattice [4] and its generalizations. For the sake of simplicity we shall write 
down the fermionic representation only for g n (U), where n is an even number. It is 

g n (U) = n!(n|e^M E ^'>° ^/«/0>^(i- VW|0> (57) 

where 

+oo 
— oo 

and where H(t) and H(t*) are given by (40) and sets of variables t := (h,t2, . . . ) and t* := (t-i, f_2, • • • 
are related to the set U := (Uq, Ui, . . . ) via 

oo 

Ui = 2 J2 im tm (58) 

where two semi-infinite sets of parameters t m are called higher times of the coupled charged BKP 
hierarchy. Let us figure out the equidistant energy levels case: 

-Ut = -2h (59) 

which is described by all t m vanish except t\. 

Notice that (58) is not one-to-one correspondence, given sets t and t uniquely define the variables U 
but not vice versa. 
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One need to make some comments on (57). We obtain it by developing the exponential term in the 
middle of the vacuum expectation value into the Taylor series. Then, as one can see, only ^n-th term 
contributes. After one uses certain tricks invented in [40] he gets (57). 

In the same way, using results of [40] we obtain that e~ F « is also a tau function. We shall not write 
precise formulas for this case. Thus, for special values of electric charges, q = 1, \/2 (and also for the 
case q = 2 which we omit in the present paper) we relate Coulomb particles to tau functions. 

2.6 Random turn walk in an external potential. Decay of the step func- 
tion 

Let us describe a certain model. Consider the random turn walk related to the weighted graph 
given by 

A(U)=A + (U)+A_(U), (60) 

and where _ 

A+(U) =Y, eUi - 1 ~ Ui fifi-i> A ~( U ) = Z^'^/i-i/i, (61) 

We use the same letters U = {U} as in the previous subsection in order to pay attention to 
a certain similarity of the answers for probabilities to find particles in given configurations ob- 
tained in rather different problems: particles in the thermal equilibrium and particles subjected 
to the chosen random motion. 

According to the previous consideration our model is the following. There is the infinite 
graph, which is the one-dimensional lattice Z (which we view as a set of integer points on a 
vertical line), with all neighboring sites linked by pairs of oppositely directed arrows. Each site 
may be either empty, or filled by at most one particle. The 'bottom' limit of this graph is fully 
packed by particles, while the 'top' limit is free of particles. At each time step, in the 'middle' 
of the graph we have some configuration of particles which are subject to random motion. 

In each unit time interval a particle chosen at random hops to the neighboring position 
provided it is free. It may hop either upward, or downward with different probabilities. The 
rate of the hop of a particle from the site % — 1 to the site % is given by 

up : r(i) = e' 17 ^- 1 (62) 

while the rate of the hops from the site % to the site i — 1 is given by 

down: r *(i) =e~ u ^ +u < (63) 

which is the number inverse to (62). 

One may say that we consider random turn walking particles in an external field, which 
affects the rates of the hops. Notice that this model should be considered as a generalization of 
the so-called random turn model which was introduced in [15]. (See e.g., [8], [18], for other 
approaches.) 

A simplifying auxiliary picture may be the following. Each particle moves in the potential 
field given by U where Ui is referred as the potential of a site i. The rate of a hop depends on 
the difference of the potentials between final and target sites as given by (62) and (63). 

The alternative picture is obtained via the Young diagram of a Maya diagram. The process 
under consideration is a random changing of the shape of a Young diagram, viewed as random 
gluing and evaporating of boxes to it in such a way that only one box is added (or removed) 
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at each time step. Each box has a 'gluing energy' U R — C/ R _i, given by the distance, R, of the 
box to the main diagonal of the Young diagram. (This distance is called a content of a box of 
a Young diagram [42]). An example is given below, where a Maya diagram of level and the 
related Young diagram are drawn (the Young diagram is drawn by bold lines). The distances 
(the "contents") of boxes marked by x on the figure are R = —1 for the lower box and R = 2 
for the upper box (these numbers coincide with the height coordinates of black balls on the first 
figure). For boxes marked by *, as going up and right we respectively obtain R = — 2, 0, 3. 




1 . Random turn walk 2. Random adding/removing a box to a Young diagram 

of particles related to the up/downward hops of particles on Maya 

on a Maya diagram. diagram. At unit time instant either a box has to be added 

at any of vacant places marked by star, or a box marked 

by x has to be removed. 

Figure 7: Two realizations of random turn walk 

Let us note that the one directed process of only adding of boxes (obtained from the general 
case by the imposing of condition that all Ui <C is well studied, then the problem is 

equivalent to the enumeration of the so-called standard Young tableau [42], see below. 

Let us evaluate transition probabilities. 

For the sake of simplicity we take Dirac level n to be zero and denote U\(n) by U\ where 
we use notation (43). 

First of all let us note that the weight of the transition from an initial configuration |A') to 
the final configuration | A) in T steps is given by the formula 

Wa^a(t) =e^'-^N A , v (T) (64) 

where Na,a'( t ) i s me number of ways to come from the position A' to the position A in T steps. 
This follows from the choice of the weight of a single hop in form (62)-(63) and from the 
definition of the weight. 

For the fermions in the state of the thermal equilibrium we have the partition function 

Z = Y,e~ Ux (65) 

A 

Each Z~ 1 e~ Ux yields the probability to find fermions in a state |A). 

In our case of random turn walk, the hard core particles which start from a given initial 
configuration |A') are described via 

Zx>{T) = Y t e Ux '- Ux "\x>W (66) 

A 
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(One may call N AjA /(t) kinematic entropy factor.) Then 

Pa'-a(t) = ^(Tj-^^-^N^^T) (67) 

yields the probability to find the hard core particles in a state |A) as a result of T-step random 
turn walk of an initial configuration |A') in an external potential {Ui, i E Z}. 

In particular, in case the random turn walk starts with the vacuum configuration, then 

Z (T) = J]e-^N A)0 (T) (68) 

A 

where N A;0 (t) is the number of ways to gain the configuration |A) in T steps if the initial 
configuration is the vacuum one. 

The probability to find the hard-core particles (subject to random turn walk in a potential 
{Ui, i G Z}) in a configuration |A) is 

P ^ A (T) = Z (T)- 1 e -^N Ai0 (T) (69) 

yields the probability to find the hard core particles in a state |A) as a result of T-step random 
turn walk of an initial configuration |A'). 
As we shall show below 

N A , (T) = 2^-^-. SA (t oo ) (70) 

\ 2 r 

where s A (too) is the Schur function evaluated in a special point, too=( 1,0,0,. . . ), see Appendix 
A.3 (177). 

Let us obtain (68)-(70) via the fermionic representation. Notice that we take 

[A_(U),A+(U)] = 1 (71) 

which simplifies calculations. 

We are interested in the evaluation of 

W X ^ X (T) = r A , A ,(T) = (X\(A + (U) + A4U)y\X'} (72) 

which gives the weight of the random process where an initial configuration of hard core parti- 
cles | A') comes to the final configuration (A| in T steps. 
By (71) we have 

e z(A + +A-) = e 4 e zA +e zA_ (73) 

At the present paper we shall study the simplest case: A' = 0. Then, our model describes 
a decay of the step configuration as a function of the discrete time T. Let us take into account 
that A-(U)\0) = 0. 

By (73), by Appendix A. 5 and also by (177) we obtain 

(A|e^+ +A -)|0) = e^W • e-^too) (74) 

Now let us decompose both sides in Taylor series in z to find the right hand side of (72) 
where A' = 0. 
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For time duration T = 2m + |A| we obtain 

Wo-a(t) = t! ( JL) e-^ SA ( too ) = t!2^^ • e^s^) (75) 

which, for U = 0, gives the number of ways to get a configuration |A) starting from vacuum 
configuration |0) in T-steps (for T — |A| = 2m is even and non-negative). The index m counts 
how many times particles, which are involved in the random process, move downward. The 
left hand side vanishes if T — |A| is odd. 

Also, notice that the length of partition A can not exceed the duration of motion T. The 
equality £(X) = T is related to the motion (or, the same, to the "path" as was explained in the 
subsection 2.3) where a single (the upmost on the vacuum Maya diagram) particle moves, and 
each hop is upward one. Thus in this very case the distance it has pass (which is equal to the 
first part of the partition A, or, the same, equal to the length of the first row of Young diagram) 
is equal to T. 

Let us consider different cases described by relation (75). 

(1) For A = the Boltzmann factor e u ° = 1 and the weight is equal to the number of ways 
to return back 

WWt) = N 0l0 (T) = (t - 1)!! = e ^ T+ - (76) 

where the last relation describes the large T limit. 

(2) The case T = |A| corresponds to the non-stop forward motion: all jumps of the par- 
ticles are in upward direction, and described by Wo->a(|A|) = (A|A+'|0). If one turns to the 
description of the motion via Young diagrams, one may see that the number 

N Ai0 (|A|) = lAIMt^) =: d(X) (77) 

describes the number of ways to create a Young diagram of given shape A by gluing at random 
box by box in such a way that each time we have a Young diagram, see fig. 7. (In other words it 
is the number of the so-called standard tableau of the shape A, see Appendix A. 1 and for details 
see [42]). Thus we have 

WW|A|) = e-^|A|! SA (too) = e-^d(A), (78) 

(3) The case £(X) = 1 means that the final state is a single particle configuration (however, 
in the middle of the process more particles may be involved, thus, the problem is different from 
the random walk of a single particle restricted to the half-line, which in T — > oo limit yields 
Brownian motion on the half-line). Thus, in this case A = (Ai) where Ai is equal to the shift of 
the upmost particle from its position in the initial vacuum configuration. Suppose the hopping 
rate 

g-E/i+Ui-i = = r (79) 

does not depend on % (which was figured out in (59), r = exp(— 2^)). As a result of sim- 
ple evaluation one shows that in T — > oo limit the dominant term in the subset of i{X) = 1 
configurations A is related to the partition A = (Ai) given by 

x r l-r 2 1 / 2 (l-r 2 ) 2 \ / _ 3 x 
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The weight for one-particles configurations in large T limit near A is given by the formula which 
resembles formula for the Brownian motion 



W _ +A ~ t! exp 



(A-ry/T)' 
2v/t 



As we see the variance is given by T?. At last we note, that in case the rate (62) depends 
on site, then in a wide class of rates in large T limit A x may be evaluated as the solution of 
Ai = r(Ai)v / T. For instance, for Gauss potential, C/j = \i 2 , one obtains Ai ~ log T. 
(4) The case where T — |A| is large enough. Via Stirling's formula we obtain 



W^a(t) ^ v^(T)s e -I L/-L e^e-^s^), |A| « T (80) 




Remark 2.4. Let us notice that the sum in the right hand side of (68) may be considered as a certain 
correlation function of a discrete version of the orthogonal ensemble of random matrices. To see it 

we use the same trick as in [54]. First, according to (177) we write 

SA (too) = Hi^lZpl j h i = X i -i + T>0, i = l,...,T (81) 

Then we use the fact that we summarize a symmetric in the variables hi function multiplied by the right 
hand side of (81). In this case we can use permutations of hi, . . . , h T to replace the sum over the cones 
hi > • • • > h T > by sum over all hi > getting a factor 1 over t!. At last we obtain 

z (T) = c £ r (i + ^±1 - \ ± hY nfn) tj \ hi _ hjl (82) 

h!,---,h T >0 V i=l J i,j=l 

/ll,-,/l T >0 V 1=1 / 

= E r A.o(T) (84) 

A 

_ T 2 +T 

where c = 2 2 and where 

V^(^) = -U hi + U T -i - logT(h t + 1) + l -^h t (85) 

and E(h, q),h = (hi, . . . , /i T ) is the electrostatic energy of T charges q located on the one-dimensional 
lattice and having coordinates hi,...,h T : 

T 12 T 

E(h, q) = Y,( U hi ~ U T -i) - ^-hi - q 2 logj] \hi -hj\+ q 2 logT(hi + 1) (86) 

i=l i<j 

The number of terms in sums (82) and (83) is finite thanks to the Gamma function insertion. However 
in T — > 00 limit, the sum ranges over all positive integers K L . 

Notice that E(h, q) vanishes on the vacuum configuration A = 0. 

Terms in the right hand side of (86) are to be interpreted as follows: the probability rates given by 
U are related to the external electric potential (which may depend on a site coordinate), gives rise to 
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the first term. The third term is the Coulomb interaction of particles with the unit charge. The last term 
has a meaning of an external electric field which provides the vanishing of the Coulomb energy of the 
vacuum configuration of the particles. Two last terms appeared due to the Schur function in the right 
hand side of (75) and descended from the hard core interaction of the particles. The (discrete) time of 
the random turn walk is equated with the total number of Coulomb particles, or, the same, with the 'size' 
of orthogonal matrix whose eigenvalues are presented by the positive integers hi, . . . , h T . 

The t\,o(t) gives the weight of the random process where the vacuum initial configuration of hard 
core particles decays to a final configuration given by coordinates hi , . . . , h T in T steps. 

Now let us turn to the problem of finding large T limit. In a usual way, we variate positions 
of particles in order to obtain the dominant term in sum (82). 

First let us remind that the length, say R, of the partition we are looking for does not exceed 
T. Introducing the density of particles a < 1 (which the number of particles for one site) as a 
function of the variable hi — \ — % + R (the origin of Maya diagram is related to h = R) in a 
way 

POO 

/ a(h)dh = R (87) 
Jo 

we come to the equation for a which defines the dominant configuration of particles at time 
T >> 1 in the continues limit 



, r(h-T) n [°° a(x)dx 1, / R 2 R f°° . , J \ n 
log v : + J h^x~ + 2 g [ T + Y~2~J xa ( x ) dx )=° (88) 



h 



where P J stands for the principal value. 

Let us note that this equation may have no solutions. 

This equation one can solve in case the hopping rate (62) does not depend on site: r(i) = r, 
see (79). With the help of the formula 



2R \ - ^arcsin - l) 2x 

2 ze L dh = i g f± ( 8 9) 

x — h R 



(which may be extracted from [39], our proof the reader will find in the Appendix) we obtain 
the following solution to (88) 

a(h) = \-\ arc sin ^ - 1^ , he [0, 2R] (90) 

and 

a{h) = 1, h < 0; a(h) = 0, h> 2R, (91) 

where the length R of the partition is found by substitution of (90)-(91) into (88) and is equal 
to 



fl=2 VrT7^ (92) 

see fig. 8 which describes the decay of initial vacuum Maya diagram of level in large T limit. 

Let us consider the partition A = A(t) and it's Young diagram related to the asymptotic 
density function <j(h). Then solution (90)-(92) correspond to a Young diagram symmetrical 
under the reflection with respect to the main diagonal (this follows from the symmetry between 
particles and holes which exists in case of constant rate). The length of the partition £ (A(t) ), it's 
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area |A(t)| and the number of boxes on the main diagonal k (A(t)) are given by the following 
formulae 

l(\(T))=R = 2^0 (93) 

"2-R E>2 d D 2 t , 

(94) 



|A(T)|=/ M ft )dk - T + - = T + - = _ + 



1 + r" 



*w t » = I)/t^ (95) 

We obtain (95) using the fact that the last number is equal to the number of particles which 
passed the origin (the origin is related to h = R) after duration T. This number is obtained by 
evaluating the following integral: 

fc(A(T)) = / a(h)dh = - (96) 
Jr tt 

The derivative rffc M T )) — ^(t(1 + r~ 2 ))~^ yields current of particles through the origin related 
to asymptotic configuration. 

The number of downward steps m(r) of particles (see (75)) related to the asymptotic con- 
figuration is 

™w = z -t^ = ^75)-iVt^ <97) 

These formulae give simple answers for particular cases: (a) completely symmetric simple 
exclusion process, r — 1, (b) locking potential r — > resulting to A(t) = (c) totally asym- 
metric exclusion process (d-TASEP) where particles hop only upward, r — > oo. Two last cases 
were considered in the list of examples which follows (75). 

Let us write down a formula for the Schur function evaluated for the asymptotic partition 
A = A(t) 

Sx (t ) = e R2 (f° °W(u-loguR)du+± f f &(u)d(u')\og(uR-u' R)dudu')+- _ £ -{R 2 logR+R 2 0(l) ^gg^ 

— logT-h- 



e 2(l + r-^) 



where by dots we denote minor terms in T — > oo limit and where a(u) := a{uR). Together 
with (70) and (97) it yields the number of ways to get the asymptotic configuration in T steps 
via random turn walk in large T limit as 

Na ( t),o(t) = 2m(T ^ (T) , ^( T) (too) = T!e-^ losT -^-) losT+ - = e^+- (99) 

Notice that in the large T limit the leading term of log Na( T ),o( t ) does not depend on the rate 
r, though the asymptotic configuration A(t) depends. In the large T limit the leading term of 
log Na( T ),o(t) is the same for the locking case r — > where answer is given by (76) and in the 
d-TASEP case r — > oo where answer is given by (77). 

It will be interesting to compare considered numbers (length, weight, current,...) evaluated 
on the asymptotic configuration with their average values. 

Let us mark that the shape of the asymptotic Young diagram given by (90)-(91) coincides 
with the shape found by Kerov and Vershik [38] in the problem of study of limiting shape of 
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Figure 8: Decay of step function for the case of a constant hopping rate in T — > oo limit 



random partitions with fixed weight | Aj = AT distributed according to the so-called Plancherel 
measure. We obtain the identification if we put N = R 2 . As soon as we arrive at (75) this 
fact is not so striking because the Plancherel measure is equal to N\s\(too) 2 . Thus, it is similar 
to the variational problem for d-TASEP where (75) is replaced by (78). Then, the variational 
problem for d-ASEP may be split into two parts: given weight of partition and taking it into 
account via Lagrangian multiplier to variate the shape of the partition and obtain (90)-(91) (this 
is quite similar to the d-TASEP problem where the weight is fixed as |A| = T), then to variate 
with respect to the weight. This might be the other way to get formulae (90)-(92). 

The problem of evaluation of correlation functions for various versions of random turn walk 
is studied in the forthcoming paper [29]. 

2.7 Weights which depend on mutual configuration: coupling with Coulomb- 
type weight 

As we have seen the case U = may be interpreted as a free random turn walk of the non- 
intersecting particles on the lattice. Here "free" means that each particle hops either upward 
or downward with the same probability (provided, the target site is empty). Then, as we have 
seen, in the answer for the weight of transition between initial vacuum configuration to a con- 
figuration given by (integer) coordinates h, hard-core interaction gives rise to the factor A(h). 
Then (83) can be interpreted as a partition function for the Coulomb gas of particles (each one 
has the unit charge) living on the lattice. 

For non-vanishing U the weights of hops also depend on the location of the hop. In 
Coulomb gas picture Ui is interpreted as an electric potential at a point i. 

Now let us consider a model of random motion where the weight of each step is additionally 
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dependent on the particle configuration as follows. Consider the 'coupling' of the model studied 
in the previous subsection 2.6 to the modified Coulomb interaction (48) where particles posses 
an arbitrary charge q. Namely, suppose that the weight of the hop is additionally multiplied by 
a factor which depends on the configuration of the particles according to formula 



hi _ j-j- hi + 1 



\ 



k i + 1 h 



hj 



(100) 



/ 



for the hop of i-th particle from the site hi to the site hi + 1 (provided this site is empty), and 
according to 

i . . . v' 



hi 



V 



h n 



J=l,2,. 



hi — 1 — hi 
hi — hj 



(101) 



/ 



for the hop of i-th particle from the site hi to the site hi — 1 (provided this site is empty). 
Then instead of (75) we obtain 



W ^ x (T;q) = T! 



1 



2 m m! 



|A|-T 



e^A(too) 1+9 = t!2^ 



1 



(102) 



Normalization function is a sum of the weights over all transitions of the vacuum configura- 
tion of particles which are described by the non-intersecting random turn walk in the Coulomb 
potential of time duration T. After the corresponding changing of the weights of steps ac- 
cording to the appearance of the Coulomb interaction with a charge q, we get the resulting 
normalization function as 



Z(T,q) :=Y,M*oo)f (X\(A_ + A + Y 



AeP 



hi 



y I 

,6u r c-^ + 1) 



2 



(s) 2 n^'in*.-*. 

^ ' i=l i,j 



(103) 



(104) 



In (104) his (non-ordered) set of positive numbers, hi,..., h T , \X(h)\ := h\-\ \-h T + ^j(T- 

1), and 

(105) 



f3 = l + q 2 
V(h) = U h + f3\ogh\ 



(106) 



Thus, by choosing the rates of hops of the random turn walk we arrive to effective electric 
charges model in thermal equilibrium with particles forced by 

(1) external field given by U 

(2) hard-core interaction between particles 

(3) Coulomb interaction between particles with coupling constant q 

Note that if q 7^ we have no anymore need to mention that our particles are hard-core 
ones, since the Coulomb repulsion force does not allow to occupy the same place by more than 
one particle. 
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One can see that the special cases, q = 0,1, y/3, result in links with discrete analogues of 
orthogonal {(3 = 1), unitary (J3 = 2) and symplectic {(3 = 4) matrix ensembles [43] which as it 
is known may be related to the integrable hierarchies. 

Then, for constant rate the asymptotic configuration will be the same (90)-(91), see fig. (8), 
where now 

RlW = Vttf* (107) 

As we see the size of the asymptotic domain shrinks because of the adding of the modified 
Coulomb interaction which tries to confine particles near the vacuum configuration, see Remark 
2.1. For asymptotic configuration X(q, t) we straightforwardly obtain 

1 1 

1 / 4t \ 2(i+ 9 ^) t 1 / 4t \ i+T 1 1 / 4t \ 2(i+«? 2 ) 

fc(A(g,T)) = - — — - , m(g,T) = --- 



7r^l + r -2y -v.,-, 2 g ^ 1 + r - 2 y 4\l + r- 2 

which gives answers for the asymptotic partition length, weight, for the number of particles 
which passed the origin and for the number of backward steps. For q 2 > in the large T limit 
and we also obtain 

1 / 4t \ 
log WW = - 2TIT? y [Y+r^ ) logT + ' ' ' 

T 

logN A( , )T)i0 (T) = -logT + --- 

At last let us mark that the introduction of the free parameter q is similar to the introduction 
of a parameter s in [37], where a notion of an entropy of stochastic dynamical systems was 
studied. We hope to consider this problem in future. 

2.8 Determinantal formulae: Wick theorem and Gessel-Viemiot formu- 
lae 

Consider A+(U) of (61). Then 

±(\>\A + (U) N \\) = (A'|e A+(C/) |A>, N = |A'| - |A| (108) 
The form of right hand side allows to apply the Wick theorem: 



U m+h x -U m + hl 



(\'\e A +M\\) = det ((m\f m+h ,e A +Mf m+hi \m)) = det ^ , , 

\ J / i,j=l,...,m \ {II j — Hi) I , 

\ J / i,j=l,...,m 

(109) 

where 

hi = \i-i + m, h'i = X'i - i + m, m > £(X') > £(X) (110) 

where we use 

(0\fi + Ne A+ Wfi\0) = ^-f Ui - Ui+N (HI) 
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The right hand side also allows to consider tau functions 

r(t, U, t) = (0|e H(t) e A+(t/) e^ ( * ) |0> (112) 

as a generating function for the transition weights (\'\A + (U) N \X). 
Remark 2.5. Let us mark the formula 

l{\<) 

(\'\e A +^\\) = s y/x (t 00 )l[e- Uh ^ (113) 

i=i 

where = (1,0,0,...) and 

Xi = hi-i + m, \\ = h\-i + m, m = max(*(A),£(A')) = *(A') (114) 



The left hand side of (108) is interpreted as a iV-step random turn process, described above. 
Let us consider an example. 

Gessel-Viennot formula for binomial determinants 

Consider two sets of positive integers a k > ■ ■ ■ > a x > and b k > ■ ■ ■ > bi > 0, where 
di > % — 1, . . . , k. Following [20] let us study the following determinant 




(115) 



called binomial determinant. Let us identify sets of and a fe _j with sets of hi and h\ related 
by (110) to partitions A and A' that is 

A, = b k -i + i-k, A • = a fc _j + i- k, i = 1, . . . , k (116) 

Take 

(l fc l fc \ f 1 k 1 k \ 

iE^^E 1 *'- [y}k-^[jJ2y^2^ y2, '"i (117) 
i=i i=i / V i=i i=i / 

It is known [42] that Schur functions s\(t), where t = [x] k , vanish on partitions whose 
length exceed k: £(X) > k. 

Then, as a specification of (109) and by using formulas (31), we get that the following tau 
function 

r(k, [x] k , [y] k ) = (k\e H ^e A e R ^\k) (118) 
= (^|e H(Wfc) |A',A;)(A,A;|e ;?(Mfc) |A;)(A , ,A;|e A |A,A;) 

i(\),i(\')<k 

= E sy([x] k )s x ([y] k )det( ai Qk ), (119) 

i(\),i(\')<k 

where 

A = Y,ififi-u (120) 

is a generating function for the binomial determinants. This A is related to the motion along 
the following graph 
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There are (i+l)(i+2)(i+3) ways to 
come from a point i to the point i+3 



6 i-1 



Figure 9: The numbers to the left are the weights of the arrows 

Remark 2.6. A tau function similar to (118) was used in [27] for the study of the so-called two matrix 
model. 

The graph related to (120) is sketched out above. The weight of an arrow may be considered 
as the number of identical arrows of unit weight. Thus, the product of weights of successive 
arrows connecting two vertices yields the number of ways connecting these vertices. 

The number of ways connecting points k and i on the Maya diagram is counted by 4 = 

(OljU^'/ilO). Therefore, 

(01/^10) = (= J-^m^m =) = Q (121) 
From (1 13) we have the relation obtained in [20]: 

det(Y£*Y) =^AVA(too) (122) 

where is the so-called skew Schur function (170). 

Gessel and Viennot suggested to interpret the binomial determinant as the number of non- 
intersecting ways connecting points dj with points hi, where % = 1, . . . , k (in our approach 
connecting a pair of Maya diagrams) as it is depicted at the right hand side figure [20]. Each 
step is either down or left directed: 

Notice that if we introduce the space-time axis (x, t) as it is figured on the picture Fig. 3 
we observe non-intersecting one-dimensional random walk of k vicious walkers, moving with 
speed ±1, changing at random the direction at each time instant At — \. The particles start at 
t — their motion from positions a 1: . . . , a k . Later a particle by particle meets a Tight' signal 
spreading with the speed — 1 which starts motion in the origin at t = 0. A particle which started 
at dj meets the signal at point |6j. 

Considering random path to be an interface one can relate the right hand figure to the ran- 
domly growing area below the path. By stars we mark the last portion of area (ice covering of 
a step scala). Then the parameter k in the tau function (118) may be viewed as discrete time. 

These problems will be considered in the forthcoming paper. 
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1 . The transition weight is given by 
the number of paths to the next figure 




2. Paths connecting a pair of Maya diagrams 
may be viewed as the interface of growing area 



Figure 10: Counting connecting paths 



If we replace e A by e A ■ ■ ■ e A = e TA in the expression for tau function (1 18) we obtain T-ple 
mapping of partitions described in the next section. In certain sense the model considered in 
this section interpolates between random turn and vicious walkers models, which we consider 
in the next section. 

3 Action of GL (oo) group elements on Fock vectors and ran- 
dom processes 

Below we shall construct certain versions of stochastic motion on the space of partitions, which 
may be also called random slices ("rocks") and random layering. Namely, one may view 
stochastic dynamics on partitions as models of random slices, which may be visualized as a set 
of partitions placed one above another according to its time evolution. Each three-dimensional 
figure of height T obtained in this way present a path (13) and may be compared with a rock 
inscribed into a corner between two walls. In this section we present only descriptions of the 
models in terms of fermions without study of these models. 



3.1 Determinantal processes 

Now consider the case where each a" -1 in (6) has a special form 

w-i = gjj-i ;= e E i ,*«?r 1 /i/ fc> ( T ) = g(r) := g^- 1 ■ --g 1 ' (123) 

These exponentials are known to be elements of GL(oo) group, see [14], [32], where the im- 
portance of GL(oo) for constructing of integrable equations together with their solutions was 
established. 
Then 

r T (x, y) = (Oljfe 1 ) • • • f(x?) g(T) f( Vl ) ■ ■ ■ f(y N )\0) A^Ad/)" 1 

is known [14], [32] to be a TL tau function. In our case we want tau functions to be dependent 
on the axillary discrete variable T via (123). 
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In what follows we shall assume the following simplifying condition 

a^'- 1} = 0, t,k<0, j = 1,2,... (124) 
Let us remind definitions of the Schur function 

s„(x) = det (x h a b ) A(x)-\ hi = Ui-i + N 

and basis Fock vectors 

W,N) = f hN --- f hl \0), (v, N\ = (0\f hl • • • fh N 
Denoting the matrix elements of <?(t): 

T X , v (T) = {\,N\g(T)\v,N), (125) 
we may rewrite tau function in form 

M x >v) = ^2r x ,u(T)sx(x)s l/ (y) 

If we have (124), then Wick's theorem has especially simple form: 

</^'V' J "V" 1) > =detg"~ 1 (h^,ht 1) ), g"'\a,b) = (0\f a g™ ; ./,()} (126) 
Remark 3.1. Thanks to the Wick theorem, we write (12) as 

r T (x,y) = ({h®}) A(x)A(y) (127) 

(i) 

where the sum ranges over T sets of all admissible configurations h a ■ , j = 0, . . . , T 

P T m ({h®}) = det (S^j det (y^ T) ) jj det (V +1 ^ +1 \ fc?>)) (128) 



Now, 



p; 



r f/ Wl^ - ^ (^L>) _ ^ ({^}) 

M ^ |/la} '" r T (x,y) "Ev^T),^) (129) 

almost coincides with the formula (1.10) of [23] (under some specification of functions tp and 4> in this 
paper, and under certain notational replacements: say, of (129) are x£* of [23] and t t (x, y) of (129) 
is Zu-m of [23]), where it describes the so-called reduced probability density . The only difference 

(i) 

between formula (129) and formula (1.10) of [23] is that, in our case, variables h a take discrete values 
as they are related to the parts of the intermediate partitions, while the main example of [23] is the 
multimatrix model where these variables coincide with eigenvalues of random matrices (then, a serves 
for the number of an eigenvalue, and j serves for a number of a matrix in the chain of coupled matrices). 

More generally (129) describes a discrete analogue of the so-called multi-level determinantal en- 
sembles where the "levels" are numbered by j in (128). 

Any o?'-? -1 in form (123) gives rise to a multi-level determinantal ensemble [13]. 
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In our approach to random processes the tau function r T (x, y) plays the role of the generat- 
ing function for t\ jV (t). The normalization function is 



A 

In the present paper we are interested not in the ratio P]^ ({ha^}^ mentioned in the remark 
(see (129)) but in: 

WT) = T -l§ (130, 

yielding the normalized weight (which may be interpreted as the probability in case all weights 
are positive) to arrive to a configuration A in T steps if an initial configuration is v. 

Let us note that the most natural way to consider (128) is to use the so-called multi- 
component fermions and multi-component integrable hierarchies. This will be done below. 

3.2 Multi-component fermions 

For certain problems it is suitable to re- write (128) with the help of multi-component fermions, 

\fU) /Mr = \fU) f(*)i = n \fU) f«i — A . .A 

(see Appendix A. 2 for details). The multicomponent fermions were used in [14], [32] to con- 
struct hierarchies of multicomponent integrable equations. 

We shall use the following notations. Given set of partitions, j — 0, . . . , T, we intro- 
duce a basis Fock vector 

\^\n^;...-^\n^) := IT^ • • • n«|0, . . . , 0) 

where 

n (j) = f% ■ ■ ■ f%, =v®-i + n«, nW = €(!/?>) 



The dual Fock vector we denote by , ; ■ ■ ■ ; 



n 



(T)| 



(A<°>, i<°>; . . . ; AW, J< T W<°>; . . . ; = * A(0) ,„<■,> • • • <W 

Now, we replace (6) by 

\u, N; 0, 0; . . . ; 0, 0) -> o 1 ' ]//, AT; 0, 0; . . . ; 0, 0) — > ► o 1 ' 1 " 1 • • • o 1 ' ]//, iV; 0, 0; . . . ; 0, 0) 

(131) 

where 

N = £(u {0) ) (132) 

and where each operator o j+1 ' j is an intertwining operator between different Fock spaces 
and : 

Then (131) describes the same process on partitions as (6) conditioned by (123) and by 
(124) provided that coincides with g j+1 ' j (k, i) of (126). 
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Denoting 

G(t) = G T,T_1 • ■ ■ G 1 ' (134) 

we re- write tau function (127) as 

r T (x, y) = (N\Y(0, 0\fW( Xl ) • • • f^(x N )G(T)p\y N ) ■ ■ ■ f {1 \yi)\0, . . . , 0) A^A^ 1 

(135) 

= ^r A ,,(T)s A (x)s,(y) 

where r A)i ,(T) of (125) is reexpressed as 

r AiV (T) = (0, 0; . . . ; 0, 0; A, iV|G( T )|z/, TV; 0, 0; . . . ; 0, 0) (136) 

Tau function (135) is a specification of the multi-component tau function introduced in [14] (see 
also [34]), and almost coincides with the fermionic representation of the milti-matrix model 
presented in [31]. 

Regardless to the choice of representation for t A i ,(t) which is obtained either from the 
one-component fermionic expectation value (125), or from the multi-component one, (136), 
the formulae (130) yields the transition probability to come to a configuration A if the initial 
configuration is v. 

Example: N vicious walkers (discrete Brownian motion of N hard core particles) 
Let us apply it to the model of N vicious walkers introduced by M.Fisher in [15] (actually 
he mainly considered the so-called lock step model where all walkers start on odd (or even) 
numbered coordinates on the lattice. At each tick of the clock each walker moves either up 
or downward subject to no two walkers occupying the same site at the same time. In various 
contexts this model was studied in [6], [17], [18], [63]. 
In this case for (133) we take 

oo 

G"- 1 = exp£ {e-^^ + e-^^ fEft 1] ) (137) 

i=0 

The number j plays the role of physical time, the parameter T coincides with the duration 
of the process. Application of operator G-" -1 to a Fock vector (=Maya diagram) describes the 
unit-time transition of N particles, each of which hops to the nearest site conditioned that in 
both initial and final configurations there are no sites occupied by more than one particle. In 
(137) U plays the role of an external field: a particle located at a site % hops upward with the 
rate e u i+^ u i and downward with the rate e u i-i^ u i 

Then (131) describes the motion of hard-core particles which start with a configuration v. 
As it follows from (184) the particles are constrained to move only along the positive half-line. 
Particles are not influenced by this restriction in case their initial configuration (given by v) is 
far enough from the origin. 

The exponent in (137) should be compared with expression (61) which was used earlier to 
get the random turn motion. 

Vicious walkers on the circle 

We need to replace (137) by 

GiJ-i = exp^-i, g^-i = e-uP+rt-VfUft-i) + e -v^^ f U)fU-D (138) 
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On correlation functions. An advantage of the multicomponent approach is the simple 
representation for the correlation functions in form of determinants. 

Given set of partitions p^\ j — 1, . . . , T — 1, whose lengths, i(p^), do not exceed N, one 
may address the following questions: 

(a) what is the weight of the process where all paths (13) are constrained by the condition 
that each number pf\ i = 1, . . . ,£(p^) < N is a part of j-th intermediate partition p^ 
belonging to the path. 

Let us consider 'characteristic operators' 

^ 3 ^= X(P U) ) = ■ ■ ■ f$f$ ■ ■ ■ 039) 



where nf' = pf' — i + iV '. Introduce 

o x ( T ) = t,t-i x (t-i) t-i,t-2 . . . 2,i x (i) i,o (140) 

Then, the weight we are looking for is the following correlation function 

v f (i) (t-i)x _ (0,0;...;0 > 0;A,JV|o x (T)KJV;0,0;...;0,0) _ det (jfilf^Mf^^ 

(141) 

where hi = \ — i + N and h\ — vi — % + N (i — 1, . . . , N), and where we use the notation 
|0) := |0, 0; . . . ; 0, 0). The last equality is due to the Wick theorem. 

(b) what is the weight of the process where all paths (13) constrained by the condition that 
no one part of an j-th intermediate partition p^ belonging to the path coincides with any of 
parts of the p^\ where j = 1, . . . , T — 1. 

Let us consider 'characteristic operators' 

X U) '-=xip {i) )=f^-f^f^-f^ (142) 



where h\ 3> = pf' — i + N. Introduce 

0x (T) = T,T-l£(T-l) T-l,T-2 . . . o 2,l^(l) o l,0 (143) 

The weight of this process is the following correlation function 

K ^ - (»^MK(T)M,..,0) _ ^t((0|il! ) o,(T)^|0)) 

where h { = X { — % + iV and h\ — ^ — % + N (i — 1, . . . , N). The last equality is due to the 
Wick theorem. 

Equivalent representation is given by 

(1) _ (T _ 1} _ (0, 0; p^\n^; pW , #; A, N\o(t)\u, N; p^\n^ p^n^O, 0) 

)- ^ 

(145) 
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where n^- 1 ) = £(p^) 

This form allow to write it via the Wick theorem as 



(1) _( T _i)n _ det((0|/ fei o(t)AJO)) 



K^ x (p (L >,...,p 



By an analogy with the so-called spectral correlation functions in random matrices, we call 
the weights K v ^\(p^\ . . . , p^ 1-1 **) and K v ^\{p^\ . . . , p^^) correlation functions. 



3.3 Random layering. Chains of Darboux transformations 

Below we consider simple application of certain GL^ transformations, sometimes called Dar- 
boux transformations. For U = below and for the case of growing partitions these examples 
are mainly reformulation of different known random models considered in literature. 
(I) Random lay 

Let us consider four basic examples of GL(oo) operators. They are basic in the sense that 
their non-vanishing matrix elements o a ,a' are related to pairs of Young diagrams different by 
adding/removing layers (strips) to Young diagrams: vertical or horizontal ones. Thus, we pick 
up four different types of o, say, = exp % — 1,2, 3, 4. Let us consider these four cases 
separately. 

(1) The exponential of 

A^(x,U):=-A^(-x,U) = - J2 J2 e ~ Ul+Ul ~ k (-^ k f^ < 146 > 

fc=l,2,... i£Z 

yields laying of a given Young diagram by a vertical strip, see fig. 11, as it is given by 

o* 1 ' (*,£/) | A'>=5>> oV,(x,U) 

A 

where 

o^,(x,U) := (A|o«|A'> = e^'-^IAHA'l 

is non-vanishing only if the difference between A and A' is the so-called vertical strip [42]. 
Example: 

Now turn to the related Maya diagrams which describes ID configuration of the related 
lattice gas. Let us compare it with the random turn walk, see Fig.7 in the subsection "Decay 
of the step function". In the random turn walk model in each time instant only one particle 
hops either one step upward or one step downward to the nearest neighboring site. The Fig.l 1 
above describes the one step of a group of randomly chosen particles upward to the nearest 
neighboring sites which occurs at one time instant, this hop being conditioned that after the 
hop each site is occupied by no more than one particle. One may say that Fig. lib describes 
'trains' of particles which hop one step upward. The potential U gives rise to the gluing rate in 
2D picture and the hopping rate in ID picture which site-depending. The weight to absorb n 
boxes is proportional to x n . 

(2) The exponential of 

A W(x, U):=- J2 T, e ~ Ui+Ui+k (- x ^f^ ( 147 > 

fc=l,2,... iei 
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Fig. 1 la. Condensation of a lay onto a Youngdiagram 
A'. Random 'rain' of boxes from rightward. These boxes 
may be fitted only to admissible places. 



x 



Fig. 1 lb. Examples of diagrams of A where o\[, is nontrivial. One layer Tone strip) 
is added at random to the Young diagram A', coxes marked by x depict the 
difference between A' and A, called the vertical strip 



Figure 1 1 : A layering by a vertical strip 

yields getting away of a given Young diagram a vertical strip, that is 

oV(x,U)\\')=J2W o^(x,U) 

A 

where 

o™,(x,U) := (A|o( 2 )|A') = e^'-^'l-|A| 

is non- vanishing only if the difference between A' and A is the vertical strip [42]. 

Presenting this process as a random motion of ID particles on Maya diagram we see that it 
describes a hop downward to the nearest site of randomly chosen group of particles (they may 
form 'trains') at each time instant. 

(3) The exponential of 

a^\x,u)= J2 e ~ Ui+Ui - kxk f^-* ( 148 > 

k=l,2,... ieZ 

yields laying of a given Young diagram by a horizontal strip, as it is given by 

o^\x,U)\\') = J2 I A) o®,(x,U) 

A 
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where the transition weight 

o^,(x,U) := (A|o( 3 )|A'> = e^-^^IAI-IVl 

is non-vanishing only if the diagram A includes the diagram A', and the difference between 
these diagrams is the so-called horizontal strip [42]. 
Example is given by fig. 12: 



j Fig. 12a. Condensation of a lay onto a Young diagram 

I I A'. Random 'rain' of boxes from downward. These boxes 

may be fitted only to admissible places, which means 
that new figure should be a Young diagram again 



x 









X 




X 


X 




X 









Fig. 12b. Examples of diagrams of A where o\[, is nontrivial. One layer (one strip) 
is added at random to the Young diagram A'. Boxes marked by x depict the 
difference between A' and A, called the horizontal strip 



Figure 12: A layering by a horizontal strip 

Turning to a ID configuration of the related lattice gas we see that it is suitable to inter- 
change roles of particle (pictured as black balls) and free sites (white balls). In this dual picture 
trains of holes hop one step downward at each time instant. The rate of this process depends on 
the difference of potential in initial and final positions and is proportional to x n , where n is the 
number of absorbed boxes. 

(4) The exponential of 

A^(x,U): = J2 e ~ Ui+Ui+kxk f^ < 149 ) 

fe=l,2,... ieZ 

yields the random process of getting away of a given Young diagram a horizontal strip, that is 

o^(x,U)\\i) = J2\V o[%(x,U) 

A 
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where 



of x ,{x,U) := (A|o( 4 )|A') = e^-^^IVI-IAI 



is non-vanishing only if the difference between A' and A is the horizontal strip [42]. 

These are four cases. Now one can consider a chain of transformations, each transformation 
is given by a set of corresponding U and x which define gluing rates. 

Examples of these chains were considered in the reviewing paper [65] from the point of 
view of discrete Hirota equations. 

(II) Chains of Darboux transformations. Random layering and "rocks" 

The chain of transformations of initial Fock vector is as follows: 

| A') ^e ACT1 ^|A') ^e^^-.-e^^lA') (150) 

where Oj = 1, 2, 3, 4, and where T is a discrete time. Given j = 1, 2, 3, . . . , T, each A" 3 (xj) 

is characterized by a given set of {uf\i G Z} which enters definitions (146),(147),(148) and 
(149). 

Given set {aj, j = 1, 2, 3, . . . , t} and {U-^ , j — 1, 2, 3, . . . , T, i e Z} gives rise to a chain 
of Darboux transformations. Such chains may be viewed as words of lengths T formed by 
four types of characters (a 'DNA' coding a set of random processes of adding or eliminating of 
layers). Tau function is constructed as 

r(T,t,t) = (0|e"(%(T)e 5 «|0> = 5> A , A ,(T)s A (t)s A ,(t) 

A, A' 

o(t) = e A ^^ ■ ■ -e^ 1 ^ 

then, each word is related to a certain generalized Darboux transformation of the tau func- 
tion. This tau function may be evaluated because we know all matrix elements of each Tetter' 

e A j J ( Xj \ Let us note that tau functions of this type were considered in [52], [51] as series in 
skew Schur functions which in special cases provided examples of Gelfand-Graev hypergeo- 
metric series. 

This is a certain random process on partitions (and their Young diagrams) describing ran- 
dom adding and elimination of lays to their Young diagrams, each time step gives rise to an act 
of either creation, or elimination of lays, as it is shown in figures 11-14. 

For the process of time duration T, the relative weight of a configuration A is given by the 
ratio 

where 

Za'(t)=XXa'(t) 

A 

is the sum of weights of all processes of duration T started with an initial configuration A'. 

One may call the random process of creating and eliminating of lays as model of random 
slices, or, a "rock" inscribed into a corner between two walls, which may be viewed as a set 
of Young diagrams placed one above another as sections of 3D figure according to its time 
evolution. Each three-dimensional figure of height T obtained in this way presents a path (13). 
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Examples. Choose A' = 0, and let <x, = 1 (which means that we shall consider growing 
Young diagrams) and U- = U , for all j — 1, . . . , T. 

Let us note that for U = models considered below may be connected to various models 
related to the so-called semi-standard tableau (see Appendix A.l and for details [42]), like a 
model of polynuclear growth (for instance, see [15] for a review). 

Af\x j ,U)=A< 1 \z j ,U) = - J2 e ~ Ui+Ui - k (- x i) k f^-* ( 151 > 

k=l,2,... i& 

In this model, Xj, j = 1, . . . , T and U, i > are sets of parameters which define transition 
probabilities. Then the process may be called a version of random process of growing Young 
diagram which is as follows. The initial configuration is the zero diagram. Each time step one 
lay (a set of strips) is added at random to the previous Young diagram (in a way that new figure 
is again a Young diagram). This is a Markov process. The transition weight for each step, say, 
a step number j, where the initial configuration (related to the moment j — 1) is, say, 
and the next configuration (related to the moment j) is fi^\ is equal to 

W>-D^>(1) = e- U ^ +U ^xf U) -^ (152) 

This yields a un-normalized probability to add a lay (a strip) of a weight (i.e. number of boxes) 
equal to — ^~^\ to a Young diagram. 
One can see that for this random process 

W ^ x (T)=e- u \s x ( Xl ,...,x T ), Z (T) =Y t e- Ux s x (x 1 ,...,x T ) (153) 

A 

Since Uq — and sq — 1 one can notice notice that 1/Z (t) yields the probability to create 
zero Young diagram, A = 0, i.e. the probability to create no any nonvanishing diagram). 
Let us consider these formulae and find links to a few known random topics. 

(1) First, consider the Gauss potential Ui = and chose Xi = e ch '% where c is a constant 
and h[, . . . , h' T are ordered positive integers h[ > • • ■ > h' T . Then the weight 

Wo-a(t) = — — ^-det (e-^-N) , hi = Xi - i + T, (154) 

unexpectedly coincides with the (un-normalized) transition probability for T non-intersecting 
("ordered") Brownian particles on the line whose initial and final coordinates are positive in- 
tegers h[, . . . , h' T and hi, . . . , h T . In large T limit it may be also related to the Itsykson-Zuber 
integral (relation between this integral and the Brownian motion is quite known topic, for in- 
stance see [36]). 

Then there are few limiting cases for (153), which are simple. 

(2) The first is x± — • • • — x T — x < 1. Then using (174) we obtain 

Z (T)=e-^ 2 -^f[^ £ f[e-^f[\ hl -h 3 f, (3 = 1 

i=0^ '' h 1: ...,h T >0 i=l i,j=l 

where 

V hl = U hi - hilogx 
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Thus Z (j) is the partition function of a discrete version of ensemble of random orthogonal 
matrices, which is also a certain tau function: the tau function of the Pfaff lattice introduced 
in [3], [34]. Let us note that the ensemble of random orthogonal matrices is well studied object, 
see [63], [6], [33]. For Ui = it is a version of the well-known Laguerre ensemble. 

(3) A similar result we obtain in case we take Xj = e v/ ~ T ^ + ^ -1 )^, where 0, $ are parame- 
ters. Then with the help of (176) we obtain unitary orthogonal ensemble 



z.(t) = 2 n ^ e n n i«* - r 5 . p = 1 ^ 

where 

v hi = u hi - h^V^i 

where q-factorials (q; q) n are defined in (172), and Ui may be chosen as series analogues to 
(53) 

oo 

U. = J2e imV ^t m , teZ (156) 

(Such parametrization was used in [50] for different purposes. These three cases can be also 
compared with [54]). 

At last let us write down a group of examples (found in [42]) where Z (t) (and, therefore, 
P ^o(t) = Zo(t) X ) may be explicitly evaluated : 

(a) Take all Ui = 0. Then, there is a simple formula 

T T 

Z o(t) = SA ( Xl ' • • • ' x t) = Y[(l - Xj)' 1 - XiXj)- 1 

A j=l i<j 

(b) Take all e Ui — 1, % < m and ef = 0, % > m. Then 

Z (T) = sx(xt, ...,x T ) = D m /D , D m = det (xf +m ~ l - xj" 1 ) 

A 

where sum is going over all partitions whose Young diagrams are contained in the rectangle 

T X 171. 

(c) Take all e Ui — 1, % < m and ef = 0, % > m, as before, and put Xj = e ( 2T - 2 i+ 1 )</ > . Then 

T „<t>(m+2j-l) i T 2<b(m+i+j-l) i 

zb(x) = e ^ ,2t - w . ■ ■ ■ . «*> = n n ^-.) : 1 = e ^ ei * 

A j=l i<j 

where first sum is going over all partitions whose Young diagrams are contained in the rectangle 

T X 171. 

The second sum ranges over all symmetric plane partitions whose diagram are contained 
in the box T x t x in. (Diagram of a plane partition is a 3D figure constructed from identical 
cubes, whose sections by horizontal planes are Young diagrams, see [42]). Here N(Q) is the 
number of symmetric plane partitions of weight 1 1 . Detailed explanations of the right hand 
sides of formulae in the last three examples (together with some other different examples of 
particular cases of summation formulae for J2\ s a) may be found in [42] from where they were 
borrowed by the authors. 

Let us also note that the model of growing partition via layering (151) is directly related to 
the problem of enumerating of the so-called semi- standard tableau. 
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4 Discussion 



In the present paper we want to figure out links between classical integrable systems and ran- 
dom system. These links in a natural way suggest to introduce site-dependent weights for hops 
of the particles and provide certain analogies between hard-core particles in equilibrium and 
non-equilibrium states. Along this line we introduced a version of discrete ASEP with site de- 
pending hopping rates. We show that the normalization function for probabilities of this model 
is related to classical integrable hierarchy of type B rather than to the hierarchy of the type A 
(the Toda lattice hierarchy) we started with. The normalization function plays the role similar 
to the role of partition function for statistical ensembles. We present a compact formula which 
describes the decay of the step function which also converts this problem to the evaluation 
of a certain correlation function for a discrete version of orthogonal matrix ensemble, where 
potential is defined by the site-dependent hopping rates of our model. 
The further project aims the following problems: 

(1) to study a model which is a modification of the model considered in the subsection 
2.6. Namely, we replace the operator A + + A_ by a linear combination of the so-called Vira- 
soro generators Li, L and L_i which form SL{2) algebra. This model is also quite solvable, 
however we will show that it presents different behavior 

(2) to present a description of models of stochastic motion which may be obtained from 
various integrable hierarchies (various realizations of hierarchies of type A,B,C and D and 
their multi-component versions) [30] 

(3) to study the asymptotic behavior of transition probabilities between different configura- 
tions and of certain correlation functions [29] 

(4) In particular for random turn decay of step function to evaluate the time dependence of : 
(a) the mean number of involved particles, < £(X) > (b) the mean height < |A| > (c) the mean 
flow through the origin < k(A) > (k(A) is the number of hooks, or, the same, the number of 
Frobenius coordinates 

(5) to consider d-ASEP with various boundary conditions. To study the phenomenon of 
shock waves and phase transition of the first order known for the problem with open boundaries 
[55], [58] 

(6) to study the phenomenon of thermalization which should occur for finite graphs 

(7) to understand links with other approaches to the ID non-equilibrium phenomena, in 
particular with these of papers [21], [22], [59], [57]. 

(8) to understand relations with the Bethe anzats method for the ASEP model 
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A Appendices 

A.l Partitions and Schur functions 

Partitions. Polynomial functions in many variables are parameterized by partitions. A parti- 
tion is any (finite or infinite) sequence of non-negative integers in decreasing order: 



The numbers A, are called the parts of the A. The number of the parts is the length of the A, 
denoted £(X). The sum of the parts denoted | A|, is called the weight of A. If | A| = n, we say that 
A is a partition of n. It is often convenient to extend the A/s to an infinite sequence {Aj}j eN 
where 



The zero partition, with £(X) = is denoted by 0. The set of all partitions, including 0, is 
denoted by P. 

The Young diagram of a partition is defined as the set of points (or nodes) (i, j) E Z 2 such 
that 1 < j < Aj. The Young diagram is viewed as a subset of entries in a matrix with /(A) 
rows and the Ai columns, with the nodes denoted by squares aligned adjacently to form rows 
of length Ai, A 2 , . . • , stacked so that successive rows downwards have equal or diminishing 
lengths, as in the example below 



which is the diagram of the partition A = (4, 3, 1), with weight |A| = 8 and length £((X\) = 3. 
The partition whose diagram is obtained by transposition of the diagram A with respect to the 
main diagonal is called the conjugate partition and denoted by v. 

Another notation is due to Frobenius. Suppose that the main diagonal of the diagram of A 
consists of r nodes (i, i) (1 < % < r). For 1 < % < r, let «j = Aj — i be the number of nodes 
in the ith row to the right of (i, i)„ and $ = Vi — i the number of nodes in the ith column below 
(i, i). We then have 



A = (Ai, A 2 , . . . , A r , . . . ) , Ai > A 2 > • • • > A r > • • • > 



(157) 



A i: =0 if i > £(A) 



(158) 



(159) 



«i > a 2 > ■ ■ ■ > a r > 

A > (3 2 > ■ ■ ■ > P r > 



(160) 



The Frobenius notation for the partition A is then 



A = (ai, • • • ,a r \Pi, ■■■,Pr) 



a\(3). 



(161) 
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The Frobenius notation may be viewed as a decomposition of a diagram A into hooks with 
corners situated along the main diagonal, the largest hook being (cti \P±), the next one (a 2 
and so on down to the smallest (a r \/3 r ). 

For example, the partition (4, 3, 1) consists of two hooks (3, 2) and (1, 0): 



and □□ (162) 

and in Frobenius notation is written (3, 1|2, 0). If A = (a\(3), then A* = (f3\a) 

If we insert the increasing sequence of integers 1, 2, . . . |A| into a Young diagram such that 
1 is in the first box, and the numbers are increasing to the right within each row and downward 
within each column, the result is called a (standard) Young Tableau. If we insert a sequence 
1, 2, . . . n with n < |A such that the numbers are nondecreasing to the right in each row, and 
decreasing downward in each column, this is called a semi-standard Young Tableau. 
Schur functions. 

The Schur functions ^(M) associated with the partition A may be viewed either as a sym- 
metric homogeneous polynomial in N variables (x±, . . . , xn), where iV may be any integer 
> \£(X) I or, equivalently a weighted homogeneous polynomial in the infinite sequence of vari- 
ables (ti,t 2 , . . .) defined by: 

N ^ N ^ N 

sa(N) := J>^,...,- J>* a ,...) := (ti,t 2) . •.,**,...) (163) 

a=l a=l a=l 

where each tj has weight i, having total weight |A|. Because of the homogeneity condition, S\ 
can only depend on the finite set of variables (t±, . . . , £|a|). Viewed as functions of the x a 's, 
they may be defined as the following ratio of determinants (Jacobi-Trudi formula) 

Sx([x]) = -%P* \ l<iJ<N (164) 

A{xi, . . .,x N ) 

where 

N 

A( Xl , ...,x N ) = Y[( Xi - Xj ) = det(x^) (165) 

i<j 

is the Vandermonde determinant. For N < | A | , this is equivalent to the following combinatorial 
definition: 

S x ([x]) = J2<---^n (166) 

where the sum is over all semi-standard Young Tableau of shape A and ^ is the number of 
times % appears in the Tableau. 

Representation theoretically, the significance of this is that ^([x]) is the character of the 
irreducible, rank A| tensor representation of U(N) or GL(N) whose symmetry properties 
are given by the Young diagram of A. (i.e. the irreducible representation consisting of ten- 
sors obtained by first symmetrizing all components labelled by the rows of the Young di- 
agram, and then antisymmetrizing those labelled by the columns. The x a 'a are viewed as 
eigenvalues of the U(N) or GL(N) group element or one simply chooses diagonal elements 
X = diag(xi, . . . , xn) and Sa(M) is the trace of the representation evaluated at X.) 
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Polynomial functions in many variables, like the Schur functions, are parameterized by 
partitions. 

Consider a semi-infinite set of variables t = (t 1: t 2 , t 3 , . . .). Given partition A, the Schur 
function s\(t) is defined by 

oo oo 

s A (t) = det(/i Ai _ i+j (t)) 1 <. J ^ (A) , where £V/i fc (t) = exp ^ z m t m , (167) 

k=0 m=l 

and, for k < 0, we put h k = . The h k (t) is called the elementary Schur function. 
If 

1 n 

t = t(xW) = (d^W), t 2 (xW), . . . ), t m (xW) = -Vxr, (168) 

then definitions (167) and (164) are equivalent [42]: 

s A (t(xW)) = s A (a;W). (169) 

From definition (167) it follows that s A (t(a;( n ))) = if £(X) > n. 

We use the underline in s x only to distinguish the two definitions. If an n x n matrix X has 
eigenvalues x±, . . . , x n , we may denote s x (x\, . . . , x n ) by s\(X), without underline, since in 
this paper the Schur function with uppercase argument is used only in this sense. 

Skew Schur function is defined as follows 

s x/lt = det (/lA^-i+jOi^n > ( 17 °) 

Properties of these polynomials are described in details in [42]. 
Schur functions evaluated at special points. 

We need notations: 

PV(j, -(- 777,1 

(a)x := (a) Xl (a - l) Aa • • • (a - k + l)\ k , (a) m : = — — , (171) 

(q a ; q)x ■= (q a ; ^((T 1 ; <?)a 2 • • • (q a ~ k+1 ; q)x k , (q a ; q) m := (1 - q a ) ■ ■ ■ (1 - q a+m ~ l ) , 

(172) 

where k = £(X). We set (a) = (<? a ; <z)o = 1 and (a)_ fe = (g a ; g)_ fc = for k > 0. 
We introduce the following notations [54]: 

^ = (1,0,0,0,...), (173) 
t(a,l)=(^,...), (174) 

t(oo,g) = (ti(oo,g),t 2 (oo,g),...), t m {oo,q) = 1 , m = l,2, (175) 

m(l — g" 1 ) 

1 _ 

t(a,g) = {t 1 (a,q),t 2 {a, q),...) , t m {a,q) = — — , m = l,2, ... (176) 

m(l — g m J 

Note that t(a, g) tends to t(oo, q) (resp. t(a, 1)) as a — > oo (resp. g — > 1). As for too, if / 
satisfies /(c£i, c 2 t 2 , c 3 t 3 , . . . ) = c d f(ti,t 2 , t 3 , . . . ) for some d G Z, we have h d f(t(oo, q)) — > 
/(t^) as ft := lng -> 0. Below A(/i) := ]l i < i (^ - hj). 
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Lemma 1. For a partition A = (Ai, A 2 , . . . ), let hi := n + Aj — i (1 < i < n), where n > £(X)). 
Then 

SA (too) = , (177) 

SA(t(a ' 1)} " ELA! fj r(a-t + i) ' (178) 

rii=i(9; 

Afo ft ) ™ 

Note that those quantities (177)-(180) are independent of the choice of n > £(X). We also 
mark that for integer positive a the Schur functions ( 1 78) and (180) vanish if£(X) > a. 



A.2 Fermionic Fock space, gl (oo) and GL(oo) 

The following is a summary regarding the one and two-component free fermion algebra based 
on the introductory section of [28]. The reader may refer to [32], [14] for further details. 

In the following, A denotes the complex Clifford algebra over C generated by charged free 
fermions {/,, fi}i e z, satisfying the anticommutation relations 

Ifi, fj]+ = tfi, = 0, Ifi, fj} + = Sir (181) 

where [, ] + denotes the anticommutator. 
Elements of the linear part 

W : = (® meZ C/ m ) © (® meZ C/ m ) (182) 

will be referred to as a free fermions. The fermionic free fields 

/(*):=£;/***, f(y) ■= hy'"' 1 ' 

kez kez 

may be viewed as generating functions for the fj, ffs. 

This Clifford algebra has a standard Fock space representation F and dual space F (see e.g. 
[26,28]) which contain unique vacuum states |0) and (0| respectively satisfying the properties 

/ m |0) = (m<0), / m |0) = (m>0), 

(0|/ ro = (m>0), (0|/ m = (m<0). (184) 

The Fock spaces F and F are mutually dual, with the hermitian pairing defined via the linear 
form (0| |0) on A called the vacuum expectation value. This satisfies 

(0|1|0) = 1; (0|/ m / m |0) = 1, m<0; (0|/ m / m |0) = 1, m > 0, (185) 
(0|/„|0) = (0|/ n |0) = (0|/ m / n |0) = (0|/ m / n |0) = 0; (0|/ m / n |0) = 0, m + n, .(186) 
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Wick's theorem implies that for any finite set of elements {w k G W}, we have 

(0\ Wl ---w 2n+1 \0) = 0, ^ 

(0\lU 1 ■ ■ ■ W 2n \0) = ^ S 9 na ( Q \ W a(l)W a ( 2 )\0) ■ ■ ■ (0\w a{2 n-l)W a (2n)\0). (187) 

Here a runs over permutations for which er(l) < cr(2), . . . , <r(2n — 1) < a(2n) and cr(l) < 
(j(3) < ••• < a(2n- 1). 

If {wi}i=i,...,Ar, are linear combinations of the fj's only, j G Z, and {u>j}i=i,...,7v linear 
combinations of the /j's, j G Z, then(187) implies 

(0|«;i • • -w N w N - ■ -wi\0) = det ((OlwiWjlO)) |ij=i,...,jv (188) 

Following [14], [32], for all N G Z, we also introduce the states 

(7V| := (0\C N (189) 



where 



Cat 
Cat 
Cn 



and 
where 



Cn 
Cn 



= f ...f N _ 1 if AT > (190) 

= f_ 1 ...f N HN<0 (191) 

= 1 ifiV = (192) 

\N) := C N \0) (193) 

= f N -i---fo ifiV>0 (194) 

= f N ...f_ 1 if AT < (195) 

= 1 ifJV = (196) 



The states (189) and (193) are referred to as the left and right charged vacuum vectors, respec- 
tively, with charge N. From the relations 

(0\f N - k f( Xi )\0)=x?- k , (0\f- N+ k-if(y t )\0)=y?- k , A; = 1,2,... iV, (197) 

and (188), it follows that 

(N\f(x 1 ) ■ ■ ■ f(x n )\0) = 8 ntN A N (x), NeZ, (198) 
(-N\f( yi ) ■ ■ ■ f(y n )\0) = 5 n , N A N (y), NeZ. (199) 

For free fermion generators with \x\ ^ \y\, 

(0|/(x)/( 2 /)|0> = — (200) 
x-y 

Note that the expression on the right hand side is actually defined, by (185), as the infinite series 
Yl™=o y""x~ n ~ x which converges only inside \x\ < \y\. However one can consider expression 
(200) for the whole region of x and y (when \x\ ^ \y\) in the sense of analytical continuation. 
From Wick's theorem it follows that 

(n - m\f( Xl ) ■ ■ ■ f(x n )f(Vi) ■ ■ ■ f(Vm)\0) = A »Wf m M (201) 

11 »=i.-.» {Xi — Vj) 

j — l,...,m 
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A.3 Commuting flows, r functions and Schur functions 

In the theory of integrable systems the following gl(oo) operators are important: 

+00 

H m = 771 = ±1 ' ±2 ' • • • ( 202) 

i=— 00 

These operators form Heisenberg algebra relations 

[H m , H n ] = m5 m+nfl (203) 
These oscillator algebra properties together with 

H m \n) = 0; {n\H- m = 0, m>0,nGZ (204) 

allow to refer vectors (0\e H ^ and H (t) as generalized coherent states which depend on param- 
eters t = (ti, t 2 , . . . ) and t = (ti, i 2 , ■ ■ ■ ). Here 

00 00 
if (t) = ]T H m t m , H(t) = H-Jm (205) 

m=l m=l 

We have the following fermionic representation of Schur functions (which follows, as usual 
from Wick's theorem) 

Lemma 1 [14] For a x > ■ ■ ■ > a k > 0, (3-y > ■ ■ ■ > f3 k > the next formula is valid: 

(0|e H(t) /-A-i •••/-A-i/c •••/«! |0> = (-l) ft+ - +/3fc+fc S A(t) , (206) 

where in the Frobenius notation A = (a 1: . . . , a k \(3i, . . . , 
Defining, as in eq. (61) 

A +(U) ^^e^-^M-i, (207) 

we have the following expression for the matrix elements of its exponential in terms of the 
Schur function evaluate at the special value too := (1,0...) 
Lemma 2 [52] 

(\,n\e A+ W\n) = e - u ^s x {t OQ ) (208) 

This relation follows from 

A+ (U) = e H °WH_ ie - H °V\ H (U) := ]T UJJ t - ^ UJJ, (209) 

i<0 i>0 

Indeed, by definitions of H (U) and (A, n\ one evaluates: 

(X,n\e Hoiu) = c n e- u ^ n) (X,n\ (210) 
where U\(n) is given by (43), where c n is defined by 

(n\e Ho{u) \n) = c n , c n = e - u °"~ u n-\ n > 0; c n = e u - 1+ - +u ", n < (211) 
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Then, (208) follows (210) and (206). Different detailed proof was written down in [26], where 
the notation r of (62) was used. In papers [50], [26] instead of e~ Ux ^ of the present paper the 
notation r\{n) was used. 

Skew Schur function is defined as follows 

s x/fl = det (h Xi _ H _ i+j )^.^ n , (212) 

The following is a generalization of (208) in terms of skew Schur functions 
Lemma 3 [52] 

(\,n\e A +W\v,n) = e u ^~ u ^ s^t^) (213) 

A.4 Derivation of (90) and some other formulae of subsection 2.6 

Formula (90) results from 

f 2 \ - 1 arcsin (y - 1) r , 

/ 2 E vy Ld y = \og2u, ue 0,2 (214) 

Jo u-y 

where we put u = -|. We have obtained formula (214) as follows. First, it is widely used in 
random matrix theory (the proof may be found, say, in [43] or in [17]) the relation called the 
Wigner semi-circle law, which is the first relation of 

l p f 1 v^ZZ^-^o, xe[-l,l], P f —1 — = (215) 
vr y_i x-v J-i Vl - v 2 x - v 

while the last relation is obtained from the first one via taking the derivative with respect to x 
and via integrating by parts. Now, keeping in mind (arcsin v)' = (1 — v 2 )~i and taking the 
derivative of the 1. h. s. of (214) with respect to u, using ((u — y)~ l ) u = —((« — y)~ l ) y and 
integrating by parts after a shift of the integration variable (y — 1 = v) we obtain 



U 7T J —l V 1 — V 2 U 1 



where as u — 1 e [—1, 1] the second term vanishes due to (215), while the first term (which 
originates from the lower boundary term; upper boundary term vanishes) coincides with the 
derivative of the r.h.s. of (214) with respect to u. Thus we have proven (214) up to a term, say 
C, independent of u. Let us find it. Consider 



■I 



' \ - \ arcsin (y - 1) f 2 f I 1 W d \ 

dy = p 7T — arcsm \v - 1 )\ -~r lo s \ u - y\ d y 

u-y J \2 7i J V dy) 

= log a - l -P f '"f 1 "^ 1 ^ =: log „ - C, , e [0, 2], 
vr Jo Jl-(y-l) 2 



where we integrated by parts. We should find the last term, C, which (we know) is the constant 
we are looking for. We have 

^ I f 2 log \u — y\dy 1 f% , 

C = -P . y \ y =- log(x-sin0)#, xe[-l,l], 
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which we evaluate at point x — 1 (since it does not depends on x we are free to chose any point 
of the interval [—1,1]): 



2 2 n=l n=l N ' 

f 1 dz / ,^ 9s i \ /"^ 1 — cos</> , , „, 

= -/ — ( l-* 2 -*-l )=-/ , y # = 21ogcos^ 2 =-log2 
Jo z \ J Jo sm0 2 



which completes the proof of (214). 
Now let us obtain (94). 



R R R 2 R f 2R R 2 R 

\M = J2( hi + i - R "> = J2 hi -Y + 2^ / ha(h)dh- — + - 
i=i i=i ~ " 

Then ^ ^ 

= / = = — / (1 + sm 2 (b)d(b = -R 2 

2 * R Jo Ji-(!l-i) 2 ^J-\ 4 



Now let us derive (92). We have 



1 , . . . i? 2 R 

-l g(T-|A|)=T- — -- 



which we equate to log ^ and obtain (92). 
At last let us obtain (96): 



r a(h)dh = r a(h) (-\dh=---— r = = ~ 

Jr ^ } Jr a{> \dh) " 2 7iRj R /J ^ —2 
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